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1.  INTRODUCTION 


This  handbook  provides  a  practical  means  for  the  prediction  of  pilot 
alr-to-air  visual  acquisition  performance.  The  predictions  are  based  upon  a 
mathematical  model  of  visual  acquisition  that  has  been  calibrated  using  data 
from  actual  flight  tests.  In  order  to  apply  the  model  to  a  particular 
situation,  the  user  must  understand  the  model  well  enough  to  1)  verify  that 
the  model  Is  applicable  to  the  situation  of  Interest,  2)  define  the  visual 
search  conditions  for  the  situation  of  Interest,  and  3)  correctly  Interpret 
the  results. 

This  handbook  contains  several  tables  that  Illustrate  model  results.  But 
these  pre-calculated  tables  cannot  Include  all  search  conditions  that  might  be 
of  Interest.  For  analysis  of  cases  not  provided  In  this  handbook,  one  may  use 
the  computer  programs  (SEEl  and  SEE2)  that  accompany  the  handbook.  The 
SEE  programs  are  Interactive  programs  that  provide  visual  acquisition 
predictions  according  to  user  Inputs.  They  are  written  In  the  Pascal  computer 
language  and  can  be  run  on  most  microcomputers.  SEEl  provides  a  basic 
calculation  of  acquisition  probability.  SEE2  provides  a  way  of  calculating 
aircraft  speed  limits  that  maintain  a  baseline  level  of  visual  acquisition 
performance.  SEEl  and  SEE2  share  many  basic  computational  modules. 

In  order  to  carry  out  the  calculations,  the  SEE  programs  make  certain 
assumptions  that  are  not  part  of  the  model  Itself  (e.g.,  SEE  programs  assume 
that  the  traffic  Is  on  a  collision  course  whereas  the  model  can  be  applied  to 
any  flight  path).  An  effort  has  been  made  to  keep  the  programs  reasonably 
simple  while  preserving  enough  flexibility  to  analyze  the  cases  of  greatest 
Interest  with  regard  to  air  safety. 
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2.  DETERMINING  APPLICABILITY  OF  THE  MODEL 


Before  using  the  model  to  analyze  a  particular  case.  It  should  be 
determined  that  the  model  Is  applicable  to  that  case*  This  section  discusses 
criteria  that  should  be  applied  In  order  to  establish  applicability. 

1)  Search  Should  Occur  Under  Daylight  Conditions 

The  model  applies  only  when  the  aircraft  Is  seen  In  bright  daylight 
conditions  (so  that  It  is  the  physical  structure  of  the  aircraft  that  must  be 
acquired,  not  the  aircraft  lights).  The  model  has  not  been  applied  to  the 
detection  of  strobe  lights,  contrails,  or  smoke. 

2)  A  Value  for  6  Must  be  Available 

In  order  to  apply  the  model,  a  search  effectiveness  parameter,  3,  must  be 
selected.  3  Is  the  rate  of  visual  acquisition  per  solid  angle  of  target  size 
per  second  of  search.  It  serves  as  as  a  curve-fitting  parameter  that  accounts 
for  all  the  human  performance  effects  upon  search  effectiveness.  3  can  be 
expected  to  vary  with  workload,  pilot  visual  acuity,  pilot  training,  and  the 
amount  of  effort  devoted  to  visual  search*  The  model  obtains  an  appropriate 
value  of  3  directly  from  flight  test  data*  However,  flight  test  data  Is 
available  for  only  a  limited  number  of  situations.  If  the  search  conditions 
of  Interest  do  not  seem  to  correspond  well  with  any  documented  flight  test, 
then  the  only  way  to  select  3  is  by  extrapolation.  The  selection  of  3  is 
discussed  In  Section  6.1. 

3)  Unusual  and  Non-standard  Visual  Conditions  Must  be  Excluded 

The  model  does  not  provide  any  explicit  adjustments  for  phenomena  such  as 
empty  field  myopia,  glint,  hypoxia,  monocular  viewing,  etc.  Any  such 
phenomena  present  In  the  flight  tests  will  be  reflected  In  the  value  of  3  that 
Is  produced.  However,  If  such  a  phenomena  Is  absent  In  the  flight  test  data, 
then  the  model  provides  no  defined  way  of  taking  It  Into  account. 

4)  Aircraft  Must  Approach  on  Non-accelerated  Collision  Courses 

Although  the  general  model  can  be  applied  to  arbitrary  flight  paths,  the 
SEE  programs  assume  that  the  target  aircraft  Is  on  a  collision  course.  This 
Is  the  situation  of  greatest  Interest  in  evaluating  the  reliability  of  visual 
separation.  The  calculations  are  also  restricted  to  rectilinear 
(unaccelerated)  flight  paths.  For  such  cases,  the  range  (separation  between 
aircraft)  decreases  at  a  constant  rate  to  zero.  This  establishes  a  one-to-one 
correspondence  between  tlme-to-colllslon  and  range.  Furthermore,  the  approach 
bearing  and  the  visual  area  presented  by  each  aircraft  remain  constant. 

5)  The  Bank  Angle  and  Pitch  Angle  of  the  Target  Aircraft  Must  be  Small 

This  Is  because  the  module  for  computing  the  visual  size  of  the  target 
assumes  that  the  aircraft  Is  seen  from  the  zero-elevation  plane.  Hence,  the 
upper  surface  of  the  wings  and  fuselage  Is  not  Included  In  the  calculation  of 
visual  area.  If  the  bank  angle  or  pitch  of  the  target  aircraft  Is  not  zero, 
the  visual  area  of  the  target  will  be  underestimated.  See  Section  4.2.3  for  a 
description  of  the  visual  area  calculation. 
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6)  Cases  for  which  the  Aircraft  Remain  Near  the  Visual  Resolution 
Threshold  for  Long  Periods  of  Time  Must  be  Excluded 

The  model  employs  a  sharp  cut-off  in  search  effectiveness  when  target 
size  is  below  a  specified  resolution  limit  (nominal  value:  1  minute  of  arc 
in  diameter).  If  the  size  is  slightly  below  the  resolution  cut-off  assumed  in 
the  model,  then  the  model  will  predict  that  the  target  will  never  be  acquired, 
no  matter  how  long  the  pilot  searches.  If  the  size  is  slightly  above  the 
resolution  cut-off,  then  the  model  will  predict  that,  given  enough  time,  the 
target  is  certain  to  be  acquired.  These  all-or-nothing  results  may  fail  to 
reflect  actual  experience  since  a  single  threshold  does  not  apply  to  all 
pilots  and  the  threshold  may  actually  vary  with  time.  Fortunately,  this 
limitation  is  almost  never  a  practical  concern  for  aircraft  on  collision 
courses  since  the  visual  target  Increases  steadily  in  size  and  is  normally 
well  above  the  resolution  limit  before  visual  acquisition  is  likely. 

7)  The  Model  Should  Not  Be  Used  to  Resolve  Fleld-of-Vlew  Issues 

The  SEE  program  provides  a  simple  test  for  whether  or  not  the  approach 
bearing  of  the  target  is  within  the  pilot's  f ield-of-view  (FOV).  The  FOV  test 
is  described  in  detail  in  Section  4.2.4.  It  is  Intended  only  to  provide  a 
rough  liidicatlon  of  the  geometries  for  which  cockpit  visibility  limitations 
will  prevent  visual  acquisition.  It  does  not  model  visibility  limitations  in 
elevation,  nor  does  it  consider  the  placement  of  window  pests  or  other 
possible  obstructions.  Thus,  it  will  not  detect  FOV  limitations  such  as  might 
arise  when  one  aircraft  approaches  the  other  from  far  above  or  far  below  the 
horizontal  plane.  The  existing  FOV  test  is  useful  in  providing  a  general 
correction  for  the  effects  of  cockpit  visibility  when  a  variety  of  speeds  and 
geometries  are  averaged.  If  it  is  necessary  to  apply  the  model  to  a  specific 
situation  where  cockpit  visibility  is  an  issue,  the  exact  fleld-of-view  should 
be  determined  using  other  techniques. 

8)  Conclusions  Must  Not  Be  Based  Upon  an  Extremely  Rapid  Increase  in 
Acquisition  Probability  that  Occurs  Just  Prior  to  Collision 

In  cases  where  visual  acquisition  is  very  difficult,  the  predicted 
probability  of  visual  acquisition  remains  very  low  until  a  few  seconds  before 
collision  and  then  rapidly  climbs  to  unity.  This  is  because  the  SEE 
calculations  assume  that  search  effort  is  constant  and  that  the  range  will  go 
to  zero  at  collision.  Under  these  assumptions,  the  visual  size  of  the  target 
becomes  so  large  that  visual  acquisition  is  certain  to  occur  before  zero  range 
is  reached,  even  if  only  in  the  last  tenth  of  a  second  before  collision.  In 
reality,  two  things  prevent  this.  The  first  is  that  the  range  to  the  pilot's 
eyes  may  not  be  zero,  even  at  collision.  But  more  Importantly,  visual  s  iarch 
effort  is  not  really  continuous  (as  the  model  assumes),  but  occurs 
intermittently  as  the  pilot's  visual  scan  moves  from  Inside  to  outside  the 
cockpit.  The  actual  probability  of  acquisition  will  lag  the  model  prediction 
during  periods  when  the  predicted  probability  is  rapidly  Increasing.  Thus, 
the  model  prediction  that  the  probability  of  acquisition  will  shoot  up  to 
unity  in  the  last  few  seconds  before  collision  is  an  artificial  result  for 
most  actual  flight  conditions.  This  is  not  normally  a  problem,  since  most 
analyses  will  discount  visual  acquisitions  that  occur  so  late  that  effective 
avoidance  is  unlikely  to  be  achieved. 
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3.  GENERATING  A  SEEl  ANALYSIS 


This  section  provides  the  Information  needed  to  run  the  basic  visual 
aqulsltlon  program,  SEEl.  A  description  of  the  computational  techniques 
employed  In  SEEl  is  provided  In  Section  4. 

SEEl  computes  tables  containing  the  probabilities  of  visual  acquisition, 
PACQ,  for  a  subject  aircraft  that  encounters  a  target  aircraft  under  specified 
visual  search  conditions.  The  search  conditions  are  specified  by  a  set  of  16 
Input  variables.  If  tables  for  several  sets  of  conditions  are  of  Interest, 
they  must  be  computed  as  separate  cases. 

For  some  purposes,  multiple  runs  of  SEEl  are  necessary.  For  example.  If 
the  probability  of  either  aircraft  acquiring  Is  desired,  then  SEEl  must  be  run 
twice,  with  one  aircraft  being  the  subject  In  the  first  run  and  the  target  In 
the  second  run.  SEEl  does  not  provide  averaging  over  several  tables,  nor  does 
It  provide  random  selection  of  Input  variables. 

3. 1  Row  and  Column  Variables 


The  PACQ  tables  are  created  by  varying  two  user-selected  Input  variables. 
The  variable  that  is  assigned  anew  for  each  row  is  called  the  row  variable. 

The  variable  that  Is  assigned  anew  for  each  column  Is  called  the 
column  variable.  The  other  14  variables  are  constant  for  all  entries  in  the 
table.  For  example,  SEEl  could  produce  a  table  of  PACQ  values  for  various 
values  of  subject  airspeed  and  visual  range.  For  all  values  in  this  table, 
other  variables  (such  as  the  airspeed  of  the  target  aircraft)  would  be 
constant. 

The  methods  used  to  assign  values  to  the  row  and  column  variable 
differ.  The  row  variable  is  stepped  from  a  minimum  value  to  a  maximum  value 
with  a  fixed  step  interval.  The  column  variable  is  stepped  through  fixed 
values  Input  by  the  user.  Up  to  8  values  are  permitted.  They  need  not  be 
spaced  evenly  or  be  monotonlc. 

Figure  3.1  depicts  the  basic  input/output  structure  of  SEEl.  The  sixteen 
(16)  Input  variables  that  must  be  specified  In  order  to  compute  a  single  value 
of  the  acquisition  probability  (PACQ)  are  listed  as  variables  1  through  16  in 
Table  3.1.  Each  variable  has  a  nominal  value  that  Is  employed  unless  the  user 
specifies  otherwise.  In  addition  to  the  16  Input  variables,  there  are  certain 
Internal  variables  that  are  Important  at  Intermediate  stages  of  the 
calculation.  These  variables  are  also  listed  In  Fig.  3.1  and  Table  3.1.  The 
details  of  the  computational  modules  are  discussed  In  Section  4.  Figure  3.2 
Illustrates  how  the  16  Input  variables  can  be  determined  from  a  description  of 
a  particular  encounter  scenario. 

3.2  Input  Variables 

The  following  paragraphs  provide  a  brief  description  of  each  of  the  16 
Input  variables.  A  guide  to  selecting  proper  values  of  variables  is  provided 
in  Section  6. 
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INPUT 


OUTPUT 


Figure  3. 1  Basic  input/output  structure  of  SEEl .  Sixteen  input  variables  must 
be  defined  to  compute  the  cumulative  probability  of  visual  acquisition  (PACQ). 
See  Table  3.1  for  variable  definitions.  See  Fig.  4.1  for  a  more  detailed  diagran 
of  the  modular  structure. 
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TABLE  3.1 


VARIABLES  USED  IN  PROGRAM 


NO. 

SYMBOL 

NAME 

NOMINAL 

VALUE 

DEFINITION 

1 

- 

BETO 

17000  per 
steradlan-sec 

Single-pilot  8  value  before  transition. 

2 

- 

BETl 

17000  per 
ateradian-sec 

Sxngle-pllot  8  value  after  transition. 

3 

tl 

Tl 

180  sec 

Modified  tau  time  for  8  transition 
(time  to  closest  approach). 

4 

t2 

T2 

12  see 

Modified  tau  time  fox  evaluation  of  PACQ 
(time  to  closest  approach). 

5 

D1 

D1 

0.  nml 

Modified  tau  range  parameter  for  8 
transition. 

6 

D2 

D2 

0.  nml 

Modified  tau  range  parameter  for  PACQ 
evaluation. 

7 

R 

R 

20.  nml 

Visual  range. 

6 

d 

DLIM 

1  arc-min 

Resolution  limit  of  human  eye. 

9 

Vl 

VSUBl 

180  kt 

Airspeed  of  subject  aircraft. 

10 

V2 

VTAR 

130  kt 

Airspeed  of  target  aircraft. 

11 

- 

lACT 

1 

Aircraft  type  for  target  aircraft 
(see  table). 

12 

X 

XANG 

180.  deg 

Crossing  angle  (relative  heading). 

13 

PL 

PL 

1.0 

Search  intensity  for  pilot  In  left  seat. 
(0  “  no  pilot  1  ■■  normal  search) 

14 

PR 

PR 

0.0 

Search  Intensity  for  pilot  In  right  seat^ 
(0  -  no  pilot  search,  1  ■  normal  search) 

15 

FL 

FOVL 

-120  deg 

Field-ot-vlew  cutoff  on  left  side  for 
searching  pilot  In  left  seat  (measured 

positive  clockwise  from  straight  ahead). 
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TABLE  3.1 


VARIABLES  USED  IN  PROGRAM  (CONT’D) 

NCMINAL 


NO. 

SYMBOL 

NAME 

VALUE 

DEFINITION 

16 

FR 

FOVR 

90  deg 

Field-of-vlew  cutoff  on  right  side  for 
searching  pilot  In  left  seat  (measured 
positive  clockwise  from  straight  ahead). 

17 

• 

r 

RDOT 

380  kt* 

Range  rate  of  target  aircraft. 

18 

AX 

35  sq.ft.* 

Head-on  visual  area  of  target. 

19 

h 

AY 

85  sq.ft.* 

Broadside  visual  area  of  target. 

20 

Az 

AZ 

260  sq.ft.* 

Topside  visual  area  of  target. 

21 

A 

AREA 

35  sq.ft.* 

Pinal  visual  area  of  target. 

22 

0^ 

THl 

0  deg* 

Bearing  of  target  aircraft  as  seen  from 
subject  aircraft. 

23 

®S 

TH2 

0  deg* 

Bearing  of  subject  aircraft  as  seen  from 
target  aircraft. 

24 

Ri 

RNGl 

3.78  nml* 

Range  at  which  8  transition  occurs. 

25 

R2 

RNG2 

1.27  nml* 

Range  at  which  PACQ  is  evaluated. 

26 

^0 

BEFO 

17000*  per 

Effective  value  of  6  before  transition. 

steradlan-sec 

27 

ei 

BEFl 

17000*  per 

Effective  value  of  8  after  transition. 

steradlan-sec 

28 

p 

PACQ 

- 

Probability  of  visual  acquisition. 

* 


Value  that  is  computed  when  all  16  Input  variables  have  nominal  value. 
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ENCOUNTER  DESCRIPTION 


CONCLUSIONS 


•The  subject  aircraft  is  flying  at  a  speed  of  •Vl  =  180kt 

180  knots. 

•  The  prevailing  visual  range  is  20  nmi.  •  R  =  20  nrai 

•There  is  (me  irilot  in  the  left  seat  and  (Hie  in  •PL=1.0.  PR=1.0, 
the  right  seat  Ncmnal  visual  acuity  ig)plies.  dt  .tm  =  i  arc-min 

•  The  cockpit  has  a  synunetric  fleld-of-view.  .  POVL  =  120  deg,  FOVR  =  90  deg 

The  pilot  in  the  left  seat  can  see  all  apiRoadi 

bearings  from  8  o'clock  around  to  3  o'clock  (i.e., 
ftom  120  degrees  to  the  left  of  the  nose  of  the 
aircraft  to  90  degrees  to  the  right). 

•  The  pilots  are  unaware  of  the  presence  of  .  beTO  =  17000/ster-sec  (the  value  for 

the  tar^  aircraft.  Woricload  is  low  and  pilots  are  unalerted  seiaxdi) 

expending  a  ncHmal  amount  of  effort  on  visual 
sesuch. 

•  A  target  aircraft  begins  to  i9>proach.  It  is  a  •IACT=1 
Piper  PA-28  flying  on  an  unaccelerat^  collision 

course. 

•It  is  flying  at  a  speed  of  130  knots.  •V2=130kt 

•  It  is  on  a  head-on  course  (i.e.,  its  heading  •  XANG  =  180  deg 
differs  by  180  degrees  ftom  that  of  the  subject 

aircraft). 

•The  crew  of  die  sulyect  aircraft  receives  a  •  T1  =  40  sec,  D1  =  0.3  sec 
traffic  advisi^  from  their  TCAS  unit  when  the 
traffic  is  project^  to  be  40  seconds  from  reaching  a 
range  of  03  nmi. 

•  Both  CTew  members  then  concentrate  on  .  beTI  =  130000/ster-sec  (the 

finding  the  indicated  traffic.  value  for  alerted  search) 

•  Visual  search  is  considered  to  have  failed  •  T2  =  12  sec,  D2  =  0.0  nmi 
unless  acquisition  occurs  prior  to  12  seconds  ftom 

collision 


Figure  3.2  Parameter  Selection  Example 


(1)  Input  Variable  BETO. 

SEEl  allows  the  single-pilot  value  of  model  parameter  B  to  transition 
once  from  an  early  value  to  a  final  value*  This  Is  used  to  model  the  effect 
of  receiving  a  traffic  advisory.  BETO  Is  the  early  value  of  B,  prior  to 
receipt  of  any  traffic  advisory.  See  further  description  of  the  transition 
logic  In  Section  4.2.6. 

(2)  Input  Variable  BETl 

This  Is  the  final  value  of  B  that  applies  after  the  transition  event. 

If  It  Is  set  equal  to  BETO,  then  the  timing  of  any  transition  does  not 
affect  visual  search.  See  further  descrlpton  In  Section  4.2.6. 

(3)  Input  Variable  T1 

The  B  transition  occurs  T1  seconds  prior  to  the  time  the  aircraft  will 
reach  a  range  of  Dl.  See  further  discussion  In  Section  4.2.6. 

(4)  Input  Variable  T2 

The  probability  of  visual  acquisition  (PACQ)  output  by  SEEl  Is  the 
cumulative  probability  of  acquisition  at  T2  seconds  prior  to  the  target 
aircraft  reaching  a  range  of  02  nautical  miles*  See  further  discussion  In 
Section  4.2.5. 

(5)  Input  Variable  01 

01  Is  the  range  threshold  criteria  associated  with  the  B  transition. 

It  Is  used  In  conjunction  with  T1  to  determine  the  time  of  B  transition.  See 
further  discussion  In  Section  4.2.6. 

(6)  Input  Variable  02 

02  Is  the  range  threshold  criteria  associated  with  the  point  at  which 
the  probability  of  visual  acquisition  Is  to  be  evaluated.  See  further 
discussion  In  Section  4.2.5. 

(7)  Input  Variable  R 

R  Is  the  prevailing  visual  range  (nml)  along  the  line  of  sight  to  the 
target  aircraft.  This  Is  the  range  at  which  any  target,  regardless  of  Its 
size  and  contrast,  tends  to  become  Invisible  to  the  human  eye.  It  Is  defined 
technically  as  the  range  at  which  the  atmospheric  conditions  reduce  the 
contrast  of  a  target  to  5  per  cent  of  Its  Inherent  (close  range)  contrast. 

See  Appendix  A  for  a  description  of  how  R  enters  Into  the  calculation  of 
the  acquisition  rate. 
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(8)  Input  Variable  DLIM 

DLIM  is  the  assumed  resolution  limit  (in  arc-min)  of  the  human  eye* 

SEEl  assigns  a  zero  acquisition  rate  to  targets  beyond  this  limit*  The 
nominal  value  of  DLIM  is  1*0  arc-min*  Normally,  this  variable  has  little 
Impact  upon  the  calculated  PACQ  values  since  most  of  the  opportunity  to 
acquire  a  closing  target  comes  after  the  subtended  angle  of  the  target  has 
become  substantially  larger  than  the  resolution  limit*  See  further 
discussion  in  Section  4*2*7* 

(9)  Input  Variable  VSUB 

VSUB  is  the  airspeed  of  the  subject  aircraft*  Note  that  SEEl  does  not 
Include  a  wind  model,  hence  all  variables  are  defined  in  the  alrmass 
coordinate  system* 

(10)  Input  Variable  VTAR 

VTAR  is  the  airspeed  of  the  target  aircraft* 

(11)  Input  Variable  lACT 

lACT  is  the  index  to  target  aircraft  type  as  defined  in  the  table  of 
aircraft  types  from  file  SEEAC  (see  Appendix  C)*  For  example,  IACT“1 
selects  a  PA-28  target;  IACT-2  selects  a  Boeing  727  target* 

Note  that  the  type  of  the  subject  aircraft  need  not  be  defined,  since 
the  physical  size  of  the  subject  aircraft  is  irrelevant  to  its  ability  to 
acquire  the  target*  However,  the  speed  of  the  subject  aircraft  (VSUB)  is 
needed  to  calculate  the  rate  at  which  the  target  is  closing  upon  the 
subject* 

A  new  aircraft  type  can  be  inserted  at  run  time  by  simply  typing  ”-1" 
in  place  of  a  defined  lACT  value*  SEEl  will  then  ask  for  the  name  of  the 
new  aircraft  type  and  for  its  three  principal  areas,  AX,  AY,  and  AZ  (see 
Section  4*2*3  for  discussion  of  how  the  principal  areas  are  used  in  the 
calculations)*  A  new  aircraft  type  defined  in  this  way  is  temporary,  and 
will  not  be  retained  for  use  in  subsequent  sessions*  In  order  to 
permanently  enter  a  new  aircraft  type  in  the  aircraft  list,  the  aircraft 
description  must  be  entered  in  the  disk  file  of  aircraft  types  (SEEAC)*  Any 
text  editor  that  does  not  Insert  word-processing  control  codes  into  a  text 
file  can  be  used  for  this  purpose  (see  Appendix  C). 

For  a  new  aircraft  type,  values  of  the  principal  areas  can  be  found  by 
tracing  the  aircraft  silhouettes  onto  a  sheet  of  grldded  graph  paper, 
counting  the  number  of  grid  squares  within  each  silhouette,  and  multiplying 
the  total  by  the  number  of  square  feet  in  a  single  grid  square* 
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(12)  Input  Variable  XANG 


XANG  is  the  crossing  angle  (relative  heading)  of  the  target  relative  to 
the  subject  aircraft  in  the  airmass  coordinate  system.  XANG-ISO® 
corresponds  to  a  head-on  encounter.  XANG«0°  corresponds  to  one  aircraft 
overtaking  the  other  from  the  six  o'clock  position.  See  Section  4.2.2  for 
further  discussion  of  XANG. 

(13)  Input  Variable  PL 

PL  is  equal  to  1  if  there  is  a  pilot  in  the  left  seat  of  the  subject 

aircraft.  PL  is  equal  to  0  if  there  is  no  pilot  in  this  seat. 

(14)  Input  Variable  PR 

PR  is  equal  to  1  if  there  is  a  pilot  in  the  right  seat  of  the  subject 

aircraft.  PR  is  equal  to  0  if  there  is  no  pilot  in  this  seat. 

(15)  Input  Variable  FOVL 

FOVL  is  the  bearing  of  the  left-side  visibility  cut-off  for  the  field- 
of-vlew  from  the  left  seat  of  the  subject  aircraft.  It  is  measured  in 
degrees  clockwise  from  the  nose  of  the  aircraft  (thus,  FOVL»-120‘’ 
corresponds  to  the  8  o'clock  bearing  position).  See  Section  4.2.4  for 
further  description. 

(16)  Input  Variable  FOVR 

FOVR  is  the  bearing  of  the  right-side  visibility  cut-off  for  the  field- 
of-vlew  from  the  left  seat  of  the  subject  aircraft.  It  is  measured  in 
degrees  clockwise  from  the  nose  of  the  aircraft  (thus,  FOVR«90® 
corresponds  to  the  3  o'clock  bearing  position). 

It  should  be  noted  that  the  fleld-of-vlew  parameters  (FOVL  and  FOVR)  do 
not  change  automatically  as  the  aircraft  type  is  changed.  They  must  be  set 
independently  if  different  values  are  desired  because  of  different  aircraft 
types. 


3. 3  Steps  in  Generating  a  SEEl  Analysis 

Figure  3. 3  provides  a  set-up  form  that  can  be  used  in  assembling  the 
information  required  for  a  SEEl  analysis.  A  guide  to  selection  of  values  for 
input  variables  can  be  found  in  Section  6.  Figure  3.4  provides  a  summary  of 
the  steps  Involved  in  generating  a  SEEl  analysis. 
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Figure  3.4  Srep-by-step  instructions  for  generating  a  SEEl  analysis. 


4.  SEEl  ANALYSIS  PROGRAM 

This  section  provides  a  description  of  the  computational  structure  of 
SEEl.  Although  It  Is  possible  to  run  SEEl  without  understanding  this 
structure,  the  Interpretation  of  SEEl  outputs  Is  aided  by  familiarity  with 
the  following  material. 

4. 1  Basic  Computational  Structure  of  SEEl 

SEEl  employs  a  number  of  modules  (subroutines)  that,  when  called  In  the 
proper  order,  compute  a  probability  of  visual  acquisition  (PACQ)  for  a  given 
set  of  search  conditions.  The  search  conditions  are  defined  In  terms  of  the 
set  of  16  Input  variables  shown  earlier  In  Table  3.1.  The  program  must  be 
given  all  16  Input  variables  In  order  to  compute  a  single  value  of  PACQ.  Each 
variable  has  a  nominal  value  that  is  employed  unless  the  user  Inserts  a 
non-nomlnal  value. 

The  block  In  Fig.  3.1  entitled  "SEEl  Computational  Modules"  consists 
of  seven  component  modules.  A  more  detailed  flowchart  that  shows  these 
modules  Is  provided  In  Fig.  4.1.  The  16  input  variables  are  all  shown  in 
this  figure  as  external  Inputs  to  the  modules  In  which  they  are  used.  In 
addition,  11  computed  variables  are  shown  as  outputs  of  the  modules  In  which 
they  are  computed.  These  computed  variables  are  listed  In  Table  3.1  as 
variables  number  17-27. 

One  Input  variable  Is  selected  to  be  the  row  variable.  It  is  stepped 
from  a  minimum  value  (xmln)  to  a  maximum  value  (xmax)  with  a  constant  step 
interval  (xstep).  The  column  variable  is  not  incremented  in  this  way  — 
Instead  the  user  must  type  In  the  actual  values  to  be  used.  A  maximum  of 
8  values  Is  permitted.  This  limit  Is  due  to  the  limited  amount  of  output  that 
can  be  placed  In  a  single  plot  or  In  a  single  printed  table. 

Figure  4.2  is  a  flowchart  that  depicts  the  manner  in  which  the  loops  that 
control  the  output  are  nested. 

4. 2  Description  of  SEEl  Modules 

The  computations  carried  out  In  each  of  the  seven  basic  SEEl  program 
modules  Is  described  below.  Each  module  corresponds  to  a  single  Pascal 
procedure.  A  listing  of  the  SEEl  code  for  the  modules  can  be  found  In 
Appendix  C. 


4.2.1  Module  1;  Aircraft  Type  Look-up 

Module  1  uses  the  aircraft  type  number,  lACT,  to  determine  various 
characteristics  of  the  target  aircraft.  Table  4.1  lists  the  aircraft  types 
that  are  provided  In  the  disk  file  SEEAC.AOO.  Each  aircraft  type  Is 
associated  with  an  abbreviated  name  (e.g.,  PA28)  and  three  numbers  (AX,  AY, 
and  AZ)  that  define  the  visual  area  of  the  aircraft  as  projected  upon  the 
three  principal  coordinate  planes  (see  Fig.  4.3).  A  new  aircraft  type  can 
be  created  by  assigning  table  values  from  the  console. 
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Figure  4.2  Loops  controlling  table  ouq)ut  for  SEEl. 
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TABLE  4.1 


AIRCRAFT  TYPES  USED  IN  SEE! 


lACT 

AIRCRAFT  TYPE 

NAME 

1 

Piper  PA-28 

PA28 

2 

Boeing  727 

B727 

3 

Boeing  747 

B747 

4 

McOonnell-Douglas  DCS 

DCS 

5 

Cessna  421 

C421 

6 

U-21  (King  Air) 

U21 

7 

F-18 

F18 

AX 

(sq.  ft.) 

AY 

(sq.  ft.) 

AZ 

(sq.  ft 

35 

85 

260 

400 

ISOO 

3100 

1200 

5700 

S300 

300 

1425 

2325 

81 

171 

417 

127 

267 

600 

50 

280 

540 

(Note:  New  aircraft  types  can  be  created  by  editing  the  SEEAC  file. 
See  Appendix  B  for  further  Information.) 
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,  broadside  area 


Figure  4.3  Principal  target  areas  used  to  compute  the  visual  area 
of  the  target  aircraft  as  seen  from  the  subject  aircraft. 
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4.2.2  Module  2;  Solving  the  Velocity  Vector  Triangle 


From  a  knowledge  of  the  aircraft  airspeeds  and  the  crossing  angle, 

Module  2  solves  the  velocity  vector  triangle  (Fig.  4.4)  to  determine  the 
approach  bearings  and  range  rate.  This  solu«'lon  assumes  that  the  aircraft  are 
on  a  rectilinear  collision  course  and  hence  that  the  horizontal  miss  distance 
Is  zero.  For  these  conditions,  the  values  of  the  range  rate  (RDOT)  and  the 
approach  bearings  (THS  and  THT)  are  constant. 

Approach  bearings  are  measured  positive  clockwise  from  straight  ahead  and 
are  expressed  as  a  number  between  **180  degrees  and  4-180  degrees.  The  crossing 
angle,  XANG,  Is  simply  the  relative  heading  formed  by  subtracting  the  heading 
of  the  target  aircraft  from  the  heading  of  the  subject  aircraft.  The  simple 
trigonometric  formula  that  relates  the  approach  bearings  and  the  crossing 
angle  Is 


XANG  -  180“  +  THS  -  THT 

The  approach  bearing  for  the  target  aircraft  (THT)  Is  Important  since  It 
determines  the  visual  area  that  will  be  seen  from  the  subject  aircraft  (see 
discussion  In  4.2.3).  The  approach  bearing  for  the  subject  (THS)  Is  used  In 
Module  4  to  determine  whether  or  not  the  target  Is  within  the  cockpit 
f leld-of-vlew. 


A. 2. 3  Module  3;  Computing  the  Visual  Area  of  the  Target 

The  visual  area  of  the  target  is  computed  by  projecting  the  principal 
visual  areas  AX  and  AY  onto  a  plane  that  Is  normal  to  the  llne*^f“slght 
between  the  two  aircraft.  This  plane  is  at  an  angle  THT  with  respect  to  the 
Y-Z  plane  for  the  target  aircraft.  An  approximate  correction  for  shielding  Is 
applied  by  assuming  that  the  actual  visual  area  Is  the  largest  of  the 
projections  of  AX  and  AY  plus  one-third  of  the  smallest  projection  (see 
Reference  2).  This  approximation  Is  errorless  when  the  aircraft  Is  viewed 
along  one  of  the  principal  axes.  Because  Module  3  assumes  that  the 
llne-of-slght  lies  within  the  X-Y  plane,  the  topside  area,  AZ,  Is  not  allowed 
to  contribute  to  the  target  visual  area.  This  Is  normally  a  good  assumption, 
but  can  result  In  underestimation  of  the  visual  area  If  the  pitch  angle  of  the 
aircraft  Is  large.  AZ  Is  defined  In  the  module  to  allow  for  future  extension 
of  the  software. 


4.2.4  Module  4;  Determining  the  Effective  Values  of  8 

Module  4  examines  cockpit  visibility  and  crew  complement  In  order  to 
determine  their  effect  upon  search  effectiveness.  It  does  so  by  taking  the 
nominal  single-pilot  8  values  and  summing  them  to  produce  effective  8 
values.  For  a  single  pilot  devoting  normal  effort  to  visual  search,  the 
value  of  8  before  and  after  the  transition  event  (see  4.2.6)  are  represented 
by  parameters  BETO  and  BETl.  Module  4  computes  two  corresponding  effective 
8  values,  BEFO  and  BEFl,  by  summing  the  8  values  for  each  pilot  who  Is  able  to 
contribute  to  visual  search. 
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AIRCRAFT 


Figure  4.4  The  velociw  vector  triande  is  solved  in  Module  2  in 
order  to  deteimine  the  closing  rate  and  the  bearings  of  approach. 
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The  field-of-vlew  of  the  cockpit  Is  determined  by  the  variables  FOVL  and 
FOVR  (see  Fig.  4.5).  It  Is  assumed  that  the  cockpit  visibility  Is  symmetric 
and  that  there  is  no  cut-off  in  elevation  (normally  true  for  targets  on 
collision  courses).  VThen  the  target  Is  outside  the  field-of-view,  the 
effective  B  value  Is  set  to  zero. 

By  setting  PL  and  PR,  the  user  specifies  the  number  of  pilots  that 
contribute  to  search  for  both  left  and  right  visibility  limits.  Normally  PL 
and  PR  are  either  zero  (no  pilot  In  the  seat)  or  unity  (single  pilot  in  seat 
and  searching  normally).  It  is  possible  to  use  a  value  between  0  and  1  to 
represent  a  pilot  who  is  present,  but  is  devoting  less  than  the  normal  amount 
of  time  to  visual  search.  Or  a  value  greater  than  1  can  be  used  to  represent 
search  by  an  additional  crew  member. 

For  the  nominal  input  values  given  in  Table  3.1,  it  is  assumed  that  a 
single  pilot  is  searching  from  the  left  seat  of  the  aircraft  (PL-1.0,  PR-0.0). 
For  the  nominal  values  of  crossing  angle  (XAN6-180  degrees),  the  Intruder  Is 
always  within  the  field  of  view.  Thus  for  the  nominal  Input  values,  BEFO-BETO 
and  BEFl-BETl. 


4.2.5  Module  5;  Determining  the  Range  of  Evaluation 

The  cumulative  probability  of  visual  acquisition,  PACQ,  depends  upon  the 
range,  RNG2,  at  which  PACQ  Is  to  be  evaluated.  Often  It  Is  desirable  to  set 
this  range  to  the  latest  point  at  which  sufficient  tine  exists  for  evasive 
action.  Then  PACQ  will  represent  the  probability  of  visual  acquisition  with 
sufficient  lead  time  for  avoidance. 

RNG2  is  computed  using  a  modified  tau  criterion  based  upon  two  input 
parameters  D2  and  R2: 


RNG2  -  D2  -  RD0T*T2 

From  this  equation  it  can  be  seen  that  RNG2  is  reached  when  the  target 
is  T2  seconds  from  reaching  a  range  of  D2.  With  this  formulation,  RNG2  can 
be  determined  as  either  a  constant  range,  a  constant  tlme-to-collision,  or  a 
modified  tau  criterion.  For  example,  if  T2-0,  then  the  equation  yields  a 
simple  range  criteria  for  which  PACQ  is  evaluated  at  range  D2.  If  D2-0, 
then  a  simple  time  criteria  results  In  which  PACQ  is  evaluated  when  T2 
seconds  remain  before  collision.  If  neither  T2  nor  D2  are  zero,  then  a 
modified  tau  criterion  results. 

4.2.6  Module  6;  Determining  the  Range  of  B  Transition 

In  a  particular  scenario,  the  nominal  value  of  B  for  a  single  pilot  may 
be  altered  by  an  event  such  as  the  receipt  of  a  traffic  advisory.  SEEl  allows 
the  nominal  single-pilot  value  tf  B  to  transition  from  one  value  to  another  at 
a  predefined  point  during  the  encounter.  The  initial  value  of  B  is  BETO. 

After  the  transition,  the  value  of  B  is  BETl  (see  Fig.  4.6).  The  range  at 
which  transition  occurs  is  RNGl.  It  is  computed  using  a  modified  tau 
criterion  involving  the  two  input  parameters  D1  and  Rl: 
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Figure  4.5  Cockpit  field-of-view  model  for  visual  search  from  left  and 
ri^t  seats  in  coclmit  FOVL  and  FOVR  are  the  bearing  angles  (in  degr^s 
clockwise  frx)m  12  o'clock)  at  which  the  freld-of-view  termmates  for  a  pilot 
in  the  left  seat 
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Figure  4.6  In  SEEl,  the  effective  value  of  beta  for  a  single  pilot  is  either 
zero,  BETO,  or  BETl.  The  transition  from  zero  to  BETO  occurs  at  range 
RLIM  (the  resolution  limit  of  the  eye).  The  transition  to  BETl  occurs  at 
the  range  at  which  a  traffic  alert  is  received. 
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RNGl 


D1  -  RD0T*T1 


The  transition  occurs  when  the  target  is  T1  seconds  from  reaching  a 
range  Dl.  Note  that  if  Tl-0,  then  a  simple  range  criterion  results.  If  Dl-0, 
then  a  simple  tlme-to'-colllslon  criterion  results. 

If  the  application  of  interest  does  not  require  a  transition,  the 
transition  is  easily  avoided  by  setting  BETl  equal  to  BETO.  The  settings  of 
Dl  and  T1  then  have  no  effect  upon  the  computations  (because  the  values  of  B 
are  unchanged  by  reaching  the  transition  point). 

4.2.7  Module  7:  Computing  the  Cumulative  Probability  of  Acquisition 
PACQ 

Module  7  is  the  final  module  that  actually  computes  the  probability  of 
visual  acquisition,  PACQ.  It  requires  8  Inputs  as  shown  in  Fig.  4.1. 

Module  7  assumes  that  the  acquisition  rate  drops  immediately  to  zero  when 
the  target  is  beyond  a  range,  RLIM,  corresponding  to  the  resolution  limit  of 
the  human  eye.  RLIM  is  defined  as  follows:  First,  DLIM,  AREA,  and  RDOT  are 
used  to  compute  the  range  at  which  the  solid  angle  subtended  by  the  target 
equals  that  of  a  circle  with  diameter  DLIM.  This  range  is  given  by  the 
formula 


1. 1284*SQRT(AREA)/DLIM 

where  DLIM  Is  in  radians.  Then,  If  this  range  is  greater  than  R,  the  visual 
range,  it  is  reduced  to  R.  The  net  result  is 

RLIM  -  minimum  [1. 1284*SQRT(AREA)/DLIM,  R] 

A  more  sophisticated  model  for  RLIM  could  be  constructed,  but  in  most 
analyses,  only  a  tiny  fraction  of  the  opportunity  for  visual  acquisition 
occurs  before  the  target  is  well  within  visual  range  and  well  above  the 
resolution  threshold.  Hence,  the  details  of  the  model  in  this  regard  are 
normally  of  little  consequence. 

Note  that  if  RLIM  is  less  than  RNGl,  then  8  will  already  have 
transitioned  to  BETl  before  the  target  is  close  enough  to  be  seen.  In  this 
case,  the  value  of  BETO  has  no  affect  upon  the  calculations. 

4. 3  SEEl  Analysis  Results 

Examples  of  SEEl  output  for  some  typical  cases  are  provided  in  this 
section. 

Figure  4.7  contains  SEEl  output  in  which  T2,  the  time  of  evaluation, 
is  the  row  variable.  Because  D2-0,  T2  is  the  time  to  projected  collision. 
The  analysis  shows  how  the  cumulative  probability  of  visual  acquisition 
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SEEl  VISUAL  ACQUISITION  ANALYSIS 
NOTE  : Example. 

SEEl. PAS  (Version  AOld) 

InFNam-  G:SEE1SET.A00  VFNam-  G:SEEVAR.AOO 
ACFNam-  G:SEEAC.A00  OutFNam-  G: SEEl. OUT 

INITIAL  VALUES  OF  INPUT  VARIABLES: 


BETO 

17000.0 

R 

2.00* 

PL 

1.00 

BETl 

17000.0 

DLIM 

1.00 

PR 

0.00 

T1 

180.0 

VSUB 

180.0 

FOVL 

-120.00 

T2 

6.00* 

VTAR 

130.0 

FOVR 

90.0 

D1 

0.00 

lACT 

1.00 

02 

0.00 

XANG 

180.0 

*  -  value  on  first  Iteration 


ROW  VARIABLE: 

T2  is  varied  from  6.00  to  60.0  in  steps  of  6.00 
COLUMN  VARIABLE  : 

R  assumes  values  2.00  3.00  5.00  15.00 

30.00  9999.00 

Target  aircraft  is  PA28 

AX-  35.0  sq  ft  AY-  85.0  sq  ft  AZ-  260.0  sq  ft  AREA-  35.0  sq  ft 
ROOT-  -310.0  kt  THT-  -0.0  deg  THS-  0.0  deg 
PACQ  TABLE  : 

TABLE  OF  PACQ  VALUES 


R 


TO 

2.00 

3.00 

5.00 

15.00 

30.00 

9999.00 

1  ^ 

6.0 

0.0720 

0.1082 

0.1525 

0.2204 

0.2430 

0.2684 

12.0 

0.0113 

0.0247 

0.0461 

0.0878 

0.1039 

0.1232 

18.0 

0.0020 

0.0076 

0.0184 

0.0439 

0.0549 

0.0687 

24.0 

* 

0.0025 

0.0082 

0.0235 

0.0307 

0.0401 

30.0 

* 

0.0006 

0.0036 

0.0123 

0.0167 

0.0226 

36.0 

* 

* 

0.0014 

0.0055 

0.0077 

0.0108 

42.0 

* 

* 

0.0002 

0.0011 

0.0015 

0.0022 

48.0 

* 

* 

* 

* 

* 

* 

54.0 

* 

* 

* 

* 

* 

* 

60.0 

* 

* 

* 

* 

* 

* 

avg. 

0.0085 

0.0144 

0.0230 

0.0394 

0.0458 

0.0536 

*-*-*-*  END  OF  ANALYSIS  *-*-*-* 

Figure  4.7.  SEEl  Analysis  Results.  A  PA28  target 
aircraft  flying  at  130  knots  approaches  in  a  head-on 
geometry.  The  visual  range  is  varied  from  2  to  9999  nml 
(the  latter  value  representing  unlimited  visual  range). 
The  probability  of  visual  acquisition  is  evaluated  at 
times  from  6  to  60  seconds  prior  to  collision. 
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grows  as  the  aircraft  approach  each  other.  Visual  range,  R,  Is  the  column 
variable.  Six  values  are  chosen.  The  value  of  9999  nml  was  used  to 
approximate  a  perfectly  clear  atmosphere.  All  variables  other  than  T2  and  R 
have  the  nominal  values. 

For  cases  In  which  the  target  aircraft  Is  outside  the  fleld-of-vlew  of 
all  crew  members,  SEEl  will  write  "xfov"  In  the  table  Instead  of  "0.0000". 

This  serves  to  explain  the  reason  for  zero  probabilities  that  may  appear 
unexpectedly. 

It  can  be  seen  from  the  table  that  If  the  visual  range  were  only  3  nml, 
then  the  predicted  probability  of  visual  acquisition  at  12  seconds  to 
collision  would  be  only  2.47  percent.  For  a  perfectly  clear  day,  the 
probability  would  Increase  to  only  12.32  percent.  Thus,  visual  acquisition  Is 
very  difficult  for  the  defined  encounter  conditions. 

Figure  4.8  Is  a  similar  analysis,  but  the  target  aircraft  Is  now  a  Boeing 
727  flying  at  240  knots.  The  closing  rate  Increases  by  60  knots,  but  because 
the  target  is  so  much  larger.  It  Is  easier  to  acquire.  The  probability  of 
visual  acquisition  for  a  perfectly  clear  day  is  now  63  percent  at  12  seconds 
to  collision. 

Figure  4.9  is  a  SEEl  analysis  in  which  the  row  variable  is  the  crossing 
angle  (XANG).  The  PACQ  table  shows  that  the  crossing  angle  can  have  an 
important  impact  upon  the  probability  of  visual  acquisition. 

Additional  examples  of  SEEl  analyses  are  provided  in  Figs.  4.10  through 

4.12. 
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SEEl  VISUAL  ACQUISITION  ANALYSIS 
NOTE  : Sample  SEEl  output* 

SEEl. PAS  (Version  AOld) 

InFNam-  G:SEE1SET.A00  VFNam-  G:SEEVAR.AOO 
ACFNam-  G.'SEEAC.AOO  OutFNam-  G:SEE1.0UT 


INITIAL  VALUES  OF  INPUT  VARIABLES: 


BETO  17000*0 

R 

2*00* 

PL  1*00 

BETl  17000*0 

DLIM 

1*00 

PR  0*00 

T1  180*0 

VSUB 

180*0 

FOVL  -120*00 

T2  6*00* 

VTAR 

240*0 

FOVR  90*0 

D1  0*00 

lACT 

2*00 

D2  0*00 

XANG 

180*0 

*  ■  value  on 

first  iteration 

ROW  VARIABLE: 

T2  is  varied  from 

6*00  to 

60*0  in 

steps 

of  6*  00 

COLUMN  VARIABLE  : 

R  assumes  values 

2*00 

3*00  5*00  15* 

00 

30*00  9999*00 

Target  aircraft  is  B727 

AX-  400*0  sq  ft  AY-1900. 

.0  sq  ft 

AZ-3100.0 

sq  ft 

AREA-  400*0  sq 

ROOT-  -420*0  kt  THT- 

-0*0  deg 

THS- 

0*0  deg 

PACQ  TABLE  ; 


TABLE  OF  PACQ  VALUES 


R 


T2 

6*0 

2*00 

3*00 

5*00 

15*00 

30*00 

9999*00 

0*2617 

0*4085 

0*5703 

0*7719 

0*8265 

0*8812 

12*0 

0,0285 

0*0854 

0*1894 

0*4220 

0*5170 

0*6332 

18*0 

* 

0*0191 

0*0709 

0*2445 

0*3367 

0.4660 

24*0 

* 

0*0021 

0*0286 

0*1515 

0*2312 

0.3556 

30*0 

* 

* 

0*0115 

0*0986 

0*1650 

0.2788 

36*0 

* 

* 

0*0039 

0*0665 

0*1210 

0*2226 

42*0 

* 

* 

0*0003 

0*0460 

0*0905 

0*1797 

48*0 

* 

* 

* 

0*0324 

0*0685 

0*1461 

54*0 

* 

* 

* 

0*0231 

0*0522 

0.1189 

60*0 

* 

* 

* 

0*0167 

0*0399 

0*0966 

avg* 

•0.0290 

0*0515 

0*0875 

0*1873 

0*2449 

0*3379 

end  of  analysis  *- 

Figure  4* 

8*  SEEl 

Analysis 

Results* 

A  B727 

target  aircraft 

flying  at  240  knots  approaches  in  a  head-on  geometry*  The 
visual  range  is  varied  from  2  to  9999  nml  (the  latter 
value  representing  unlimited  visual  range)*  The 
probability  of  visual  acquisition  is  evaluated  at  times 
from  6  to  60  seconds  prior  to  collision* 
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SEEl  VISUAL  ACQUISITION  ANALYSIS 
NOTE  : Sample  SEEl  output. 

SEEl. PAS  (Version  AOld) 

InFNam-  G: SEEl SET. AOO  VFNam-  G:SEEVAR.A00 

ACFNam-  G:SEEAC.A00  OutFNam-  G: SEEl. OUT 

INITIAL  VALUES  OF  INPUT  VARIABLES: 


BETO 

17000.0 

R 

2.00* 

PL 

BETl 

17000.0 

DLIM 

1.00 

PR 

T1 

180.0 

VSUB 

180.0 

FOVL 

T2 

12.0 

VTAR 

130.0 

FOVR 

D1 

0.00 

lACT 

1.00 

D2 

0.00 

XANG 

0.00* 

*  *  value  on  first  iteration 

ROW  VARIABLE: 

XANG  is  varied  from  0.00  to  180.0  in  steps  of  30.00 
COLUMN  VARIABLE  : 

R  assumes  values  2.00  3.00  5.00  15.00 

30.00  9999.00 

Target  aircraft  is  PA28 

AX-  35.0  sq  ft  AY-  85.0  sq  ft  AZ-  260.0  sq  ft 
PACQ  TABLE  : 

TABLE  OF  PACQ  VALUES 
R 


YA  MP 

2.00 

3.00 

5.00 

15.00 

30.00 

9999.00 

a/VIM  V 

0.0 

0.9726 

0.9863 

0.9934 

0.9975 

0.9982 

0.9987 

30.0 

0.8089 

0.8895 

0.9415 

0.9780 

0.9843 

0.9895 

60.0 

0.2754 

0.3934 

0.5205 

0.6816 

0.7271 

0.7739 

90.0 

0.0895 

0.1536 

0.2397 

0.3808 

0.4292 

0.4841 

120.0 

0.0338 

0.0659 

0.1148 

0.2051 

0.2389 

0.2792 

150.0 

0.0157 

0.0332 

0.0611 

0.1153 

0.1361 

0.1612 

180.0 

0.0113 

0.0247 

0.0461 

0.0878 

0.1039 

0.1232 

avg. 

0.3153 

0.3638 

0.4167 

0.4923 

0.5168 

0.5443 

*-*-*-*  END  OF  ANALYSIS  *-*-*-* 


Figure  4.9.  SEEl  Analysis  for  Varying  Crossing  Angle 
(XANG).  The  subject  aircraft  is  flying  at  180  knots. 
PA-28  target  aircraft  flying  at  130  knots  approaches  at 
angles  from  0  (subject  overtakes)  to  180  degrees 
(head-on).  The  visual  range  is  varied  from  2  to  9999 
nmi  (the  latter  value  representing  unlimited  visual 
range).  The  probability  of  visual  acquisition  is 
evaluated  at  12  seconds  prior  to  collision. 


1.00 

0.00 

120.00 

90.0 
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SEEl  VISUAL  ACQUISITION  ANALYSIS 
SEEl.PAS  Version  A03 

InFNam=  GrSEElIN.AOO  VFNam-  GrSEEVAR.AOO 

ACFNam-  GrSEEAC.AOO  OutFNam-  G: SEE, OUT 

INITIAL  VALUES  OF  INPUT  VARIABLES: 


BETO 

17000.0 

R 

2.00* 

PL 

1.00 

BETl 

17000.0 

DLIM 

1.00 

PR 

0.00 

T1 

180.0 

VSUB 

180.0 

FOVL 

-120.00 

T2 

12.0 

VTAR 

180.0 

FOVR 

90.0 

D1 

0.00 

lACT 

2.00 

D2 

0.00 

XANG 

0.00* 

*  “  Initial  value  (before  Iteration) 


ROW  VARIABLE:  XANG  is  varied  from  0.00  to  180.0  in  steps  of  30.00 
COLUMN  VARIABLE:  R  assumes  values:  2.00  3.00  5.00  15.00  30.00  9999 

Target  aircraft  is  B727 

AX-  400.0  sq  ft  Ay-1900.0  sq  ft  AZ-3100.0  sq  ft 


TABLE  OF  PACQ  VALUES 
R 


2.00 

3.00 

5.00 

15.00 

30.00 

9999.00 

XANG 

0.00 

1.0000 

1,0000 

1.0000 

1.0000 

1.0000 

1.0000 

30.00 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

60. 00 

0.  9884 

0.9993 

1.0000 

1.0000 

1.0000 

1.0000 

90.00 

0.6530 

0.8683 

0.9708 

0.9988 

0.9997 

1.0000 

120. 00 

0.2872 

0,5255 

0.7627 

0.9539 

0.9796 

0.9943 

150.00 

0.1129 

0.2523 

0.4504 

0.7459 

0.8282 

0.9060 

180.00 

0.0632 

0.1516 

0.2925 

0.5546 

0. 6471 

0.7508 

avg. 

0. 5864 

0.6853 

0.7823 

0.8933 

0.9221 

0.9501 

*-*-*-*  END  OF  ANALYSIS  *-*-*-* 


Figure  4.10.  SEEl  Analysis  for  Varying  Crossing  Angle 
(XANG).  A  B727  target  aircraft  flying  at  180  knots 
approaches  at  angles  from  0  (subject  overtakes)  to  180 
degrees  (head-on).  The  visual  range  is  varied  from  2  to 
9999  nmi  (the  latter  value  representing  unlimited  visual 
range).  The  probability  of  visual  acquisition  is 
evaluated  at  12  seconds  prior  to  collision. 
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SEEl  VISUAL  ACQUISITION  ANALYSIS 
NOTE  : Sample  analysis. 

SEEl. PAS  Version  A03 

InFNam-  GiSEElIN.AOO  VPNam-  GtSEEVAR.AOO 

ACFNam-  GrSEEAC.AOO  OutFNam-  G; SEE. OUT 


INITIAL  VALUES  OF  INPUT  VARIABLES: 


BETO 

17000.0 

R 

2.00* 

PL 

1.00 

BETl 

17000.0 

DLIM 

1.00 

PR 

0.00 

T1 

180.0 

VSUB 

130.0 

FOVL 

-120.00 

T2 

12.0 

VTAR 

130.0 

FOVR 

90.0 

D1 

0.00 

lACT 

1.00 

D2 

0.00 

XANG 

0.00* 

*  ■  Initial  value 

(before  Iteration) 

ROW  VARIABLE: 

XANG  Is 

varied  from 

0. 00  to 

180.0 

In 

steps  of 

30.00 

COLUMl 

«  VARIABLE:  R 

assumes  values  :  2.00 

3.00 

5. 

00  15.00 

30.00  9999 

Target 

aircraft 

Is  PA28 

AX«>  35.0  sq  ft 

AY-  85. 

0  sq  ft  AZ' 

-  260.0 

sq 

ft 

TABLE 

OF  PACQ 

VALUES 

2.00 

R 

3.00 

5.00 

15.00 

30.00 

9999.00 

XANG 

0.00 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

30.00 

0.9831 

0.9939 

0.  9980 

0.9996 

0.9998 

0.9999 

60.00 

0.4459 

0.5746 

0.6942 

0.8221 

0.8538 

0.8845 

90.00 

0.1525 

0.2367 

0.3392 

0.4891 

0.5367 

0.5887 

120.00 

0.0600 

0.1048 

0.1670 

0.2716 

0. 3085 

0.3509 

150.00 

0.0314 

0.0589 

0.0990 

0.1706 

0.1967 

0.2274 

180.00 

0.0234 

0.0449 

0.0767 

0.1341 

0.1552 

0.1799 

avg. 

0.3852 

0.4305 

0.4820 

0.5553 

0.5787 

0. 6045 

*-*-*-*  END  OF  ANALYSIS  *-*-*-* 


Figure  4.11.  SEEl  Analysis  for  Varying  Crossing  Angle 
(XANG).  The  subject  aircraft  Is  flying  at  130  knots.  A 
PA-28  target  aircraft  flying  at  130  knots  approaches  at 
crossing  angles  from  0  to  180  degrees  (head-on).  The 
visual  range  Is  varied  from  2  to  9999  nml  (the  latter 
value  representing  unlimited  visual  range).  The 
probability  of  visual  acquisition  Is  evaluated  at 
12  seconds  prior  to  collision. 
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SEEl  VISUAL  ACQUISITION  ANALYSIS 
S^l.PAS  Version  A03 

InFNam-  G:SEE1IN.A00  VFNam-  G:SEEVAR.AOO 

ACFNan-  G:SEEAC.A00  OutFNam-  G: SEE. OUT 


INITIAL  VALUES  OF  INPUT  VARIABLES: 


BETO  17000.0 

R 

2.00* 

PL 

1.00 

BETl  17000.0* 

DLIM 

1.00 

PR 

0.00 

T1 

180.0 

VSUB 

180.0 

FOVL 

-120.00 

T2 

12.0 

VTAR 

130.0 

FOVR 

90.0 

D1 

0.00 

lACT 

1.00 

D2 

0.00 

XANG 

90.0 

*  -  Initial  value  (before  Iteration) 

ROW  VARIABLE:  BETl  Is  varied  from  17000.00  to 

136000.0 

In  steps  of 

17000.00 

COLUMN  VARIABLE 

:  R  assumes  values 

:  2.00 

3.00  5.00  15.00 

30.00  9999 

Target  aircraft  Is  PA28 

AX-  35. 

0  sq  ft 

AY-  85.0 

sq  ft  AZ 

-  260.0  sq 

ft  AREA- 

75.7  sq  ft 

ROOT-  -222.0  kt 

THT-  -35 

.8  deg 

THS-  54. 

2  deg 

TABLE 

OP  PACQ  VALUES 

2.00 

R 

3.00 

5.00 

15.00 

30.00 

9999.00 

BETl 

17000.00 

0.0895 

0.1536 

0.2397 

0.3808 

0.4292 

0.4841 

34000.00 

0.1710 

0.2835 

0.4220 

0.6166 

0.6742 

0.7339 

51000.00 

0.2452 

0.3936 

0. 5606 

0.7626 

0.8140 

0.8627 

68000.00 

0.3127 

0.4867 

0.6659 

0.8530 

0.8938 

0.9292 

85000.00 

0.3742 

0.5655 

0.7460 

0.9090 

0.9394 

0.9635 

102000.00 

0.4302 

0.6322 

0.8069 

0.9436 

0.9654 

0.9811 

119000.00 

0.4812 

0.6887 

0.8532 

0.9651 

0.9803 

0.9903 

136000.00 

0.5277 

0.7365 

0.8884 

0.9784 

0.9887 

0.9950 

avg. 

0.3290 

0.4925 

0.6478 

0.8011 

0.8356 

0.8675 

*-*-*-*  END  OF  ANALYSIS  *-*-*-* 


Figure  4.12.  SEEl  Analysis  for  Varying  BETl.  A  PA-28 
target  aircraft  flying  at  130  knots  approaches  at  a 
crossing  angle  of  90  degrees.  The  visual  range  Is  varied 
from  2  to  9999  nml  (the  latter  value  representing 
unlimited  visual  range).  The  probability  of  visual 
acquisition  Is  evaluated  at  12  seconds  prior  to 
collision. 
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5.  SEE2  ANALYSIS  PROGRAM 


5. 1  Equivalency  Analysis 

The  visual  acquisition  model  can  be  used  to  determine  the  maximum 
aircraft  speed  at  which  a  standard  level  of  visual  acquisition  performance  can 
be  achieved.  One  use  of  such  an  analysis  might  be  to  determine  the  extent  to 
which  waivers  from  normal  speed  limits  are  justified  in  a  particular 
situation.  There  are  two  possible  approaches  to  establishing  a  maximum  speed. 
The  first  Is  to  Impose  an  absolute  requirement  on  the  achieved  probability  of 
visual  acquisition  (PACQ).  If  this  approach  Is  used,  then  the  program  SEEl 
can  be  used  to  determine  the  maximum  speed.  A  second  approach,  known  as 
equivalency  analysis,  requires  that  the  PACQ  value  be  equivalent  to  that  of  a 
reference  pair  of  aircraft  flying  In  the  same  airspace.  The  absolute  value  of 
PACQ  Is  not  fixed.  As  will  be  explained  below,  equivalency  analysis  has 
certain  advantages.  A  special  program,  SEE2,  has  been  written  to  assist  In 
equivalency  analysis. 

One  difficulty  with  the  use  of  an  absolute  PACQ  criterion  is  the  need  to 
justify  the  value  of  PACQ  that  Is  to  be  required.  A  reasonable  standard  would 
be  the  value  of  PACQ  that  results  from  compliance  with  normal  federal 
regulations  and  practices.  However,  this  value  varies  greatly  with  aircraft 
speeds,  approach  geometry,  and  visual  range.  The  PACQ  that  is  attained  under 
worst  case  legal  flight  conditions  (two  aircraft  flying  at  the  speed  limit, 
approaching  head*-on,  with  barely  legal  VFR  visual  range)  appears  to  be  too  low 
to  serve  as  a  suitable  standard.  A  more  reasonable  standard  could  be  the 
typical  PACQ  value  resulting  from  existing  regulations.  But  the  averaging 
procedure  used  to  find  the  typical  value  then  becomes  an  Issue.  Furthermore, 
the  use  of  an  absolute  standard  requires  that  the  selected  values  of  all  Input 
variables  be  accurate  in  an  absolute  sense.  Absolute  values  may  not  be  known 
for  the  traffic  environment  of  Interest.  (For  example,  PACQ  is  affected  by 
the  visual  range  that  exists.  If  the  visual  range  in  the  airspace  where  the 
waiver  will  apply  Is  uncertain,  then  the  PACQ  value  that  results  from  a  given 
airspeed  will  be  uncertain. ) 

When  equivalency  analysis  is  applied,  the  see— and~avold  performance  of 
the  subject  aircraft  against  a  typical  target  is  compared  to  the  performance 
of  a  reference  pair  of  aircraft  flying  under  Identical  flight  conditions.  The 
speed  of  the  subject  aircraft  Is  altered  until  the  see-and— avoid  performance 
for  the  pair  involving  the  subject  equals  the  see-and-avoid  performance  of  the 
reference  pair.  The  principal  output  of  the  program  is  then  the  value  of  the 
speed  required  to  produce  equivalent  confidence  In  visual  acquisition. 

The  advantage  of  equivalency  analysis  lies  in  its  Insensitivity  to  the 
exact  values  of  many  Input  variables.  This  Insensitivity  exists  because  many 
of  the  Input  variables  affect  the  reference  pair  and  the  subject  pair 
similarly.  (For  example,  each  pair  is  assumed  to  be  subjected  to  the  same 
meterologlcal  conditions,  hence  changes  in  visual  range  tend  to  have  minimal 
effect  upon  the  relative  performance  of  the  two  pairs.) 
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For  SEE2,  the  nominal  reference  pair  consists  of  two  single-engine 
general  aviation  aircraft  flying  under  visual  flight  rules.  SEE2  considers 
visual  acquisition  to  be  successful  If  either  aircraft  of  a  pair  acquires 
successfully. 

5.2  Input  Data  for  SEE2  Analysis 

Input  specification  for  a  SEE2  analysis  Is  more  Involved  thc^n  for  a  SEEl 
analysis  because  two  pairs  of  aircraft  must  be  specified.  The  Input  data  fcr 
SEE2  resides  In  an  Input  text  file.  SEE2  reads  this  file  to  Initialize  the 
analysis.  The  Input  file  can  be  altered  using  an  ordinary  text  editor. 

Figure  5. 1  Is  a  listing  of  a  typical  SEE2  input  file.  The  first  line  is 
a  file  header  that  can  be  any  string  of  alphanumeric  characters.  The  second 
Is  a  title  line  for  the  analysis.  The  third  line  defines  the  row  variable  In 
terms  of  Its  name.  Index  number,  minimum,  maximum,  and  step  size.  The  sixth 
line  defines  the  minimum  and  maximum  speeds  that  the  subject  aircraft  will  be 
allowed  to  assume.  (If  no  equivalency  can  be  found  within  these  limits,  a 
failure  to  converge  Is  declared.)  This  line  also  defines  the  precision  (In 
knots)  required  to  halt  equivalency  Iteration  (In  this  example,  0.01  knots). 
The  next  line  Is  a  heading  line.  The  next  16  lines  contain  a  table  that 
defines  the  Input  variables  for  all  four  aircraft  (two  pairs).  The  values  in 
the  first  two  columns  are  redundant,  since  the  reference  pair  must  consist  of 
two  Identical  aircraft. 

SEE2  uses  the  same  table  of  aircraft  types  as  SEEl. 

5. 3  Computational  Structure  of  SEE2 

A  flowchart  of  the  basic  structure  of  SEE2  Is  provided  In  Fig.  5.2.  A 
listing  for  SEE2  Is  provided  In  Appendix  C.  Most  of  the  basic  SEEl 
computational  modules  described  previously  In  section  4  are  employed  In  SEE2. 

A  row  variable  Is  stepped  between  defined  limits  to  generate  a  set  of  speed 
limits.  For  each  value  of  the  row  variable,  the  probability  of  acquisition  Is 
first  computed  for  the  reference  pair.  It  Is  possible  that  equivalency  cannot 
be  achieved.  In  such  a  case,  no  value  Is  generated  for  that  value  of  the  row 
variable. 

The  pair  calculation  assumes  that  PACQ  for  each  aircraft  Is  statistically 
independent.  If  pij  Is  the  probability  of  successful  acquisition  for  aircraft 
1  searching  for  aircraft  J,  then  equivalency  exists  when 

(1-P12)(1~P21)  "  (i-P34)(i-P43) 

In  this  equation,  the  left  side  represents  the  probability  of  acquisition 
failure  for  the  reference  pair.  The  right  side  represents  the  probability  of 
acquisition  failure  for  the  subject  pair. 
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DATA  FILE;  g;SEE2.S4  Version  c 
Sample  Analysis 

BETl  2  17000.  136000.  17000. 

0.  1500.  0.01  {  yll,  yul  stopltt  } 


{VARIABLE 

REF  ACl 

REF  AC2 

REF  AC  3 

SUB  AC4} 

BETO 

17000 

17000 

17000 

17000 

BETl 

17000 

17000 

17000 

17000 

T1 

180.0 

180.0 

180.0 

180.0 

T2 

12.0 

12.0 

12.0 

12,0 

D1 

0.0 

0.0 

0.0 

0.0 

D2 

0.0 

0.0 

0.0 

0.0 

R 

20.0 

20.0 

20.0 

20.0 

DLIM 

1.0 

1.0 

1.0 

1.0 

VSUB 

130.0 

130.0 

130.0 

130.0 

VTAR 

130.0 

130.0 

130.0 

130.0 

lACT 

1 

1 

1 

1 

XANG 

180.0 

180.0 

180.0 

180.0 

PL 

1.0 

1.0 

1.0 

1.0 

PR 

0.0 

0.0 

0.0 

0.0 

FOVL 

-120 

-120 

-120 

-120 

FOVR 

90 

90 

90 

90 

Figure  S. 1.  Input  Data  File  for  SEE2. 
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5.4  Examples  of  SEE2  Output 


Figures  5.3  and  5.4  provide  examples  of  SEE2  output.  In  Fig.  5.3,  the 
combined  visual  acquisition  probability  must  be  0.2676.  It  Is  seen  that 
Increasing  BETl  by  a  factor  of  7  (from  17000  to  136000)  allows  the  subject 
aircraft  to  Increase  Its  speed  from  130  kt  to  330  kt  while  maintaining  an 
equivalent  level  of  safety.  Note  that  at  130  knots,  each  aircraft  In  the  pair 
has  an  equal  probability  of  acquiring  (p>*0. 1442)  while  at  330  kt,  the  primary 
responsibility  for  visual  acquisition  must  rest  with  the  subject  (who  has 
probability  0.2418  compared  to  only  0.03402  for  the  target). 

In  Fig.  5.4,  a  similar  analysis  Is  run  for  a  reference  crossing  angle  of 
90  degrees.  Note  that  even  though  the  reference  acquisition  probability  has 
Increased  by  almost  a  factor  of  three  (from  0.2676  to  0.7622),  the  equivalent 
speeds  do  not  change  dramatically.  This  demonstrates  the  Insensitivity  of  the 
equivalency  analysis  to  the  details  of  the  environmental  conditions. 
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*-*_*-*  SEE2  VISUAL  ACQUISITION  ANALYSIS  *-*-*-* 
SEE2  VERSION:  G:8/3b/87 
INPUT  DATA  SET  NAME:  G:SEEIN.A00 
INPUT  DATA  SET  TITLE:  Sample  Analysis 
SetFNam-G:SEE2.S4  VFNam  -G:SEEVAR.AOO 

ACFNam-G: SEEAC.AOO  OutFNam-GtSEE.OUT 


INITIAL  VALUES 

OF  INPUT  VARIABLES: 

-  REFERENCE 

:  PAIR  - 

DUMMY  SUBJECT 

ACl 

AC2 

AC3  AC4 

BETO 

17000 

17000 

17000  17000 

BETl 

— 

— 

-  - 

ROW  VARIABLE 

T1 

180.0 

180.0 

180.0  180.0 

T2 

12.0 

12.0 

12.0  12.0 

D1 

0.00 

0. 00 

0. 00  0. 00 

D2 

0.00 

0.00 

0.00  0.00 

R 

20.0 

20.0 

20.0  20.0 

DLIM 

1.00 

1.00 

1.00  1.00 

VSUB 

— 

— 

-  - 

FLOATING  VARIABLE 

VTAR 

130.0 

130.0 

130.0  130.0 

I  ACT 

1 

1 

1  1 

XANG 

180.0 

180.0 

180.0  180.0 

PL 

1.0 

1.0 

1.0  1.0 

PR 

0.0 

0.0 

0.0  0.0 

FOVL 

-120 

-120 

-120  -120 

FOVR 

90 

90 

90  90 

ROW  VARIABLE:  BETl  varies  from  17000. 

00  to  136000*0  In 

steps  of  17000.00 

FLOATING  VARIABLE:  VSUB  (between  limits  of  1.00  and 

1500.00) 

AIRCRAFT  TYPES: 

ACID 

TYPE 

AX  AY 

AZ 

ACl  1 

PA28 

35.0  85. 

0  260.0 

AC2  1 

PA28 

35.0  85. 

0  260.0 

AC3  1 

PA28 

35.0  85. 

0  260.0 

AC4  1 

PA28 

35.0  85. 

0  260.0 

OUTPUT  TABLE  - 

BETl 

VSUB 

pacqref  pacq3 

pacq4 

« 

17000.000 

130.00 

0.2676  0.1442 

0.1442  » 

« 

34000.000 

175.91 

0.2676  0.0986 

0.1875  » 

« 

51000.000 

212.11 

0.2676  0.0749 

0.2083  » 

« 

68000.000 

242.26 

0.2676  0.0604 

0.2205  » 

« 

85000.000 

268.23 

0.2676  0.0506 

0.2286  » 

« 

102000.000 

291.11 

0.2676  0.0435 

0.2343  » 

« 

119000.000 

311.58 

0.2676  0.0382 

0.2385  » 

« 

136000.000 

330.12 

0.2676  0.0340 

0.2418  » 

*-*-*-*  END  OF  SEE2  ANALYSIS  *-*-*-* 

Figure  5.3.  SEE2  Output  for  Data  Set  Shown  in  Figure  5.1. 
The  allowable  speed  of  the  subject  aircraft  Increases  as 
BETl  for  the  subject  is  Increased  from  17000  to  136000 
(a  factor  of  7). 
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*-*-*-*  SEE2  VISUAL  ACQUISITION  ANALYSIS  *-*-*-* 
SEE2  VERSION:  G:8/3b/87 
INPUT  DATA  SET  NAME:  G:SEEIN.A00 
INPUT  DATA  SET  TITLE:  Sample  Analysis 
SetFNam-G:SEE2.S4  VFNam  -G:SEEVAR.A00 

ACFNam-G : SEEAC . AOO  OutFNam-G : SEE . OUT 


INITIAL  VALUES  OF  INPUT  VARIABLES: 

-  REFERENCE  PAIR  - 
ACl  AC2 


BETO 

17000 

17000 

BETl 

— 

— 

T1 

180.0 

180.0 

T2 

12.0 

12.0 

D1 

0.00 

0.00 

D2 

0.00 

0.00 

R 

20.0 

20.0 

DLIM 

1.00 

1.00 

VSUB 

— 

— 

VTAR 

130.0 

130.0 

lACT 

1 

1 

XANG 

90.0 

90.0 

PL 

1.0 

1.0 

PR 

0.0 

0.0 

FOVL 

-120 

-120 

FOVR 

90 

90 

DUMMY 

SUBJECT 

AC3 

AC4 

17000 

17000 

ROW  VARIABLE 

180.0 

180.0 

12.0 

12.0 

0.00 

0.00 

0.00 

0.00 

20.0 

20.0 

1.00 

1.00 

FLOATING  VARIABLE 

130.0 

130.0 

1 

1 

90.0 

90.0 

1.0 

1.0 

0.0 

0.0 

-120 

-120 

90 

90 

ROW  VARIABLE:  BETl  varies  from  17000.00  to  136000.0  in  steps  of  17000.00 
FLOATING  VARIABLE:  VSUB  (between  limits  of  1.00  and  1500.00) 


AIRCRAFT  TYPES: 


ACID 

TYPE 

AX  AY 

n 

ACl  1 

PA28 

35. 

0  85. 

0  260.0 

AC2  1 

PA28 

35. 

0  85. 

0  260.0 

AC3  1 

PA28 

35. 

0  85. 

0  260.0 

AC4  1 

PA28 

35. 

0  85. 

0  260.0 

OUTPUT  TABLE  -  -  - 

BETl 

VSUB 

pacqref 

pacq3 

pacq4 

« 

17000. 000 

130.00 

0. 7622 

0.5124 

0.5124  » 

« 

34000.000 

179.95 

0.7622 

0.3299 

0.6451  » 

« 

51000.000 

219.15 

0.7622 

0.2281 

0.6919  » 

« 

68000.000 

251.51 

0.7622 

0.1685 

0.7140  » 

« 

85000.000 

279.24 

0.7622 

0.1308 

0.7264  » 

« 

102000.000 

303.63 

0.7622 

0.1053 

0.7342  » 

« 

119000.000 

325.65 

0.7622 

0.0884 

0.7391  » 

« 

136000.000 

345.79 

0.7622 

0.0771 

0.7423  » 

*-*-*-*  end  OF  SEE2  ANALYSIS 

*_*-*-* 

Figure  5.4.  SEE2  Output  for  a  Reference  Crossing  Angle  of 
90  Degrees.  The  allowable  speed  of  the  subject  aircraft 
Increases  as  BETl  for  the  subject  Is  Increased  from  17000 
to  136000  (a  factor  of  7). 
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6.  SELECTION  OF  INPUT  VALUES 


This  section  provides  guidance  In  the  selection  of  Input  values*  For  some 
variables,  the  proper  value  will  be  obvious  from  the  description  of  the 
situation  of  Interest.  For  other  variables,  the  value  Is  not  obvious  without 
reference  to  flight  test  data  or  to  special  tables. 

6*1  Search  Effectiveness  Parameters:  BETO  and  BETl 


The  search  effectiveness  parameter,  ,  relates  the  rate  of  visual 
acquisition  to  the  size  of  the  target  (see  Appendix  A  for  equations).  If 
atmospheric  visibility  Is  good,  “  Is  the  acquisition  rate  that  Is  achieved  per 
steradlan  of  subtended  target  size.  For  example.  If  ^  Is  17,000/steradlan-sec 
and  the  target  subtends  1  mlcro-steradlan,  then  the  acquisition  rate  Is 
0.017/sec.  and  the  probability  of  visual  acquisition  for  each  second  of 
search  would  be  1.7  percent. 

“  Is  decreased  by  pilot  workload  and  Is  Increased  by  Increased  vigilance 
and  knowledge  of  target  location  (such  as  that  provided  by  a  traffic 
advisory).  Ideally,  °  would  be  derived  from  flight  test  data  In  which  the  key 
conditions  affecting  “  are  Identical  to  the  situation  of  Interest.  However, 
only  a  handful  of  flight  test  programs  to  characterize  visual  acquisition  have 
been  conducted.  Four  flight  tests  are  described  below. 

6.1.1  Unalerted  Search  by  General  Aviation  Pilots 

In  a  flight  test  at  M.I.T.  Lincoln  Laboratory  (reference  3),  visual 
acquisition  performance  was  measured  for  24  general  aviation  subject  pilots 
flying  a  Beech  Bonanza  on  a  cross-country  course*  Concentration  upon  visual 
search  was  avoided  by  describing  the  test  as  a  study  of  differences  In  VFR 
pilot  techniques  under  normal  flight  conditions.  Subjects  were  asked  to 
provide  periodic  workload  ratings,  answer  formal  questions,  and  to  call  out 
all  traffic  seen.  During  a  40  minute  flight,  a  second  test  aircraft  made 
three  Intercepts,  passing  500  feet  above  or  below  the  subject.  The  time  of 
each  visual  acquisition  was  recorded  and  the  range  of  each  acquisition 
determined  from  radar  tapes.  The  experimental  technique  was  Judged  successful 
In  preventing  abnormal  pilot  emphasis  upon  the  visual  search  task.  The  value 
of  "  that  resulted  was  17,000/ster-8ec* 

6.1.2  Search  by  Airline  Pilots 

In  1957,  Wayne  Howell  of  the  Civil  Aeronautics  Administration  reported 
upon  a  series  of  flight  tests  (reference  4)  In  which  subject  airline  pilots 
flew  a  DC-3  aircraft  while  being  Intercepted  by  a  second  DC-3.  A  camera  was 
set  up  to  photograph  the  eye  movements  of  the  subject.  In  one  set  of  flights, 
the  subjects  were  told  that  eye  movements  would  be  studied,  but  were  not 
Informed  that  an  Intercept  would  be  performed.  An  analysis  of  the  reported 
data  using  the  visual  acquisition  model  yields  a  **  value  of  34,000/sec*  It 
should  be  noted  that  the  emphasis  upon  eye  movements  may  have  resulted  In 
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greater  vigilance  on  the  part  of  the  aubjects  than  would  be  present  in 
non-teat  situations.  In  a  second  series  of  flights,  pilots  were  informed  that 
an  Intercept  would  be  made.  This  Information  caused  the  pilots  to  alter  their 
scanning  technique  somewhat,  but  did  not  result  in  a  clear  Improvement  in 
visual  acquisition  performance.  Thus,  the  3  value  under  both  test  conditions 
may  correspond  to  a  higher  level  of  search  effort  than  is  typical  of  normal 
flight. 


6.1.3  Alerted  Search  Using  ATARS  Traffic  Advisories 

Flight  tests  of  the  Automatic  Traffic  Advisory  and  Resolution  Service 
(ATARS)  were  conducted  at  the  M.I.T.  Lincoln  Laboratory  in  1976-77 
(reference  1).  General  aviation  pilots  flew  a  Cessna  172  aircraft  during 
planned  Intercepts.  The  traffic  advisory  display  provided  Intruder  bearing  to 
the  nearest  clock  position  (30  degree  sector).  Data  from  109  encounters  was 
analyzed  and  the  value  of  6  was  determined  to  be  90,000/ster-sec. 

6.1.4  Alerted  Search  Using  TCAS  II 

In  a  series  of  flights  at  M.I.T.  Lincoln  Laboratory  (reference  2),  visual 
acquisition  performance  was  determined  for  professional  pilots  using  the 
Traffic  Alert  Collision  Avoidance  System  (TCAS  II).  This  system  provided  a 
traffic  advisory  with  a  bearing  accuracy  of  approximately  8  degrees 
(one-sigma).  The  resulting  6  value  was  130,000/8ter-sec. 

6.1.5  Determining  6  by  Extrapolation 

When  the  conditions  of  interest  do  not  correspond  to  any  relevant  flight 
test,  the  value  of  8  must  be  obtained  by  extrapolation.  Two  general  rules  can 
be  applied.  First,  8  Increases  in  direct  proportion  to  the  amount  of  time 
devoted  to  visual  search.  Thus,  when  a  pilot  is  alerted  to  the  presence  of  a 
target,  8  can  double  even  if  the  alert  contains  no  information  on  the 
direction  from  which  the  target  is  approaching.  Second,  8  should  Increase 
Inversely  with  the  angular  size  of  the  area  that  must  be  searched.  Insofar  as 
the  search  process  consists  of  random  sweeps  throughout  the  angular  region  in 
which  the  target  may  exist,  the  time  required  to  find  the  target  is 
proportional  to  the  size  of  the  region.  More  accurate  traffic  advisory 
Information  will  Increase  8  by  reducing  the  angular  area  to  be  searched. 

6.2  Evaluation  Range  Parameters;  T2  and  D2 

T2  and  D2  determine  the  modified  tau  value  at  which  the  cumulative 
probability  of  visual  acquisition  (PACQ)  is  to  be  evaluated  (see  description 
of  Module  5  in  4.2.5).  Normally  D2  is  set  to  zero  and  T2  represents  the  lead 
time  at  which  visual  acquisition  is  required  in  order  for  visual  avoidance  to 
be  effective.  In  literature  dealing  with  visual  collision  avoidance,  a  range 
of  values  can  be  found  for  the  required  acquisition  lead  time.  The  lowest 
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value  appears  to  be  4.5  seconds  for  military  aircraft  (Reference  5).  The 
highest^^Jalue  appears  to  be  15  seconds  for  civil  transport  aircraft 
TCAS  collision  avoidance  system  (Reference  2),  Military  aircraft  are  expect 
to  Lve  a  lower  requirement  due  to  the  fact  that  1)  they  can  accelerate  more 
strongly  to  avoid  collision  and  2)  military  pilots  may  react  faster  to  a 
sighting  due  to  air  combat  training  that  has  involved  making  quick  ^neuvers 
to  avoid  objects  on  collision  courses.  The  nominal  required  lead  time  used  in 
this  handbook  are  6  s  for  military  aircraft  and  12  s  for  civil  aircraft. 

It  should  be  noted  that  when  used  as  a  requirement,  T2  is  not  intended  to 
represent  a  comfortable  lead  time,  nor  to  represent  the  smallest  time  at  which 
visual  acquisition  could  possibly  be  of  use.  Instead,  it  is  intended  to 
represent  a  point  at  which  the  probability  of  successful  visual  avoidance  has 
fallen  to  about  50  percent.  This  choice  allows  accurate  predictions  of  the 
avoidance  failure  rate  since  avoidance  failures  that  might  occur  despite 
acquisition  prior  to  T2  will  be  offset  by  avoidance  successes  that  occur 
despite  acquisition  after  T2. 


6.3  Resolution  Limit  of  Human  Eye;  DLIM 

DLIM  is  the  angular  diameter  of  the  smallest  high-contrast  circular 
object  that  can  be  detected  by  the  pilot's  eye.  From  laboratory  tests,  a 
typical  value  of  this  resolution  limit  is  1.0  arc-min.  It  is  unusual  for  any 
subject  to  achieve  much  better  than  0.5  arc-min  resolution  in  the  laboratop. 
In  an  actual  cockpit  environment,  the  effects  of  vibration,  viewing  through  a 
window,  and  so  forth  may  increase  DLIM.  In  flight  tests,  it  has  been  observed 
that  target  aircraft  are  almost  never  seen  until  they  exceed  2  arc-min.  The 
nominal  value  of  DLIM  used  in  SEEl  is  1.0  arc-min.  The  exact  value  employed 
is  seldom  of  critical  importance,  since  very  little  of  the  total  opportunity 
to  acquire  accumulates  when  the  target  aircraft  is  near  the  resolution  limit. 
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APPENDIX  A 


BASIC  EQUATIONS  OF  THE  VISUAL  ACQUISITION  MODEL 


General Formulation  of  the  Model 


This  appendix  provides  a  concise  derivation  of  the  basic  equations  of  the 
visual  acquisition  model* 

In  any  visual  search  process,  the  Instantaneous  acquisition  rate  can  be 
defined  as  the  probability  of  visual  acquisition  per  instant  of  time,  l*e*, 

P  [acq  in  At  at  time  t] 

At  ♦O  At 

The  classic  Poisson  process  models  a  situation  in  which  the  probability 
of  an  event  occurring  in  each  instant  of  time  is  constant.  In  equation  (1), 
the  acquisition  rate  can  change  with  time*  The  cumulative  probability  of 
acquisition  is  then  modeled  as  a  non-homogenous  Poisson  process.  The 
cumulative  probability  of  acquisition  by  time  t2  can  be  written 

P  [acq  by  t2]  ■  1.0  -  exp  [  -  /  X(t)  dt]  (2) 

^CP 


TABLE  A.1  NOTATION 

A  visual  area  presented  by  the  target  aircraft 

C(r)  target /background  contrast  at  range  r 

Cq  Inherent  contrast  of  target  with  background 

t  time 

r  range  (separation  between  aircraft) 

ri  range  at  which  search  begins 

r2  range  at  which  search  terminates 

R  visual  range  (point  at  which  contrast  degrades  by  95Z) 

rdot  range  rate 

&  search  effectiveness  parameter 


47 


The  above  equation  Is  a  purely  mathematical  consequence  of  the  way  in 
which  visual  acquisition  Is  described.  The  usefulness  of  the  model  hinges 
upon  our  ability  to  define  the  acquisition  rate  X  for  specified  search 
conditions.  In  examining  several  sets  of  experimental  data.  It  was  determined 
that  X  at  any  Instant  of  time  Is  proportional  to  the  solid  angle  subtended  by 
the  target,  l.e. 


A 


(3) 


It  Is  known  from  laboratory  experiments  that  It  Is  actually  the  product 
of  the  target  area  and  target  contrast  that  determine  the  detectability  of  a 
target.  Hence,  If  contrast  varies,  the  effect  can  be  modeled  by  substituting 
the  area-contrast  product  for  the  area  above: 


C(r)  A 

Co  r2 

For  vision  through  a  homogenous  atmosphere,  contrast  decreases 
exponentially  with  range  according  to  Koschmleder's  law: 


(4) 


-2.996  r 

C(r)  -  Cq  exp  [ - ] 

R 


(5) 


Thus  the  acquisition  rate  for  search  through  a  homogenous  atmosphere 
can  be  written 


A  -2.996  r 

X  «  3 - -  exp  I  - ]  (6) 

r2  R 

The  cumulative  probability  of  acquisition  by  time  t2  is  obtained  through 
Integration  of  A  according  to  equation  (2): 

t2  8A  -2.996  r 

Ptacq  by  t2]  -  1.0  -  exp  [“/  - exp  ( - )  dt]  (7) 

•  r2  R 

This  equation  Is  the  most  general  result  of  the  visual  acquisition  model 
since  it  allows  8,  A,  and  r  to  vary  with  time  In  an  arbitrary  manner. 
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Application  to  a  Nominal  Collision  Course 

For  an  aircraft  on  an  unaccelerated  collision  course.  It  can  often  be 
assumed  that  A  Is  constant  and  that  the  range  Is  decreasing  at  a  constant 
rate.  For  any  such  period  In  which  6  Is  also  constant,  the  opportunity  for 
acquisition  Is  described  by  the  Integral  of  the  size-contrast  product,  Q  : 


Q(t) 


1  -2.996  r 

- Exp  [ - ]  dt 

r2  R 


(8) 


The  ciimulatlve  probability  of  visual  acquisition  Is 

P  [acq  by  t]  -  1  -  exp  (  -  3  Q(t)I  (9) 

The  nominal  units  of  Q  are  steradlans-sec.  One  unit  can  accumulate  from  a 
target  subtending  one  thousandth  of  a  steradlan  being  In  the  field  of  view  for 
one  thousand  seconds.  In  practice,  most  targets  of  Interest  will  subtend  only 
a  few  millionths  of  a  steradlan. 

The  opportunity  that  accumulates  as  the  target  moves  from  r^  to  r2 
(r2  <  ri)  can  be  written 

A  E2(2.996r2/R]  E2(2. 996ri/R] 

|rdot|  r2  ri 

where  E2  is  the  second-order  exponential  integral 


•  Exp(-zy) 

E2(z)  -  f  - 

1  y2 

This  Integral  can  be  evaluated  by  a  series  expansion  (see  Pascal  function 
E2  in  ^pendlx  C). 

Expressions  for  Infinite  Visual  Range 

When  the  visual  range  Is  Infinite  (perfectly  clear  atmosphere),  the 
exponential  Integrals  In  (10)  have  a  value  of  unity  regardless  of  the  range  at 
which  they  are  evaluated.  The  expression  for  probability  of  acquisition 
simplifies  to 


-6A 

P[acq  by  r2]  -  1  -  exp  [  - (l/r2  “  1/ri)]  (11) 

I rdot I 
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The  probability  density  function  for  the  range  at  which  acquisition  occurs 
Is  then 


f racq(r) 


6A  -0A 

-  exp[ - (1/r  -  1/ri)]  (12) 

1 rdot I  r2  | rdot | 

0  <  r  <  ri 


0  elsewhere 

The  mean  range  of  visual  acquisition  computed  from  the  above  pdf  is 


«  (-l)nqn 

E[racqJ  -  q  ri  exp[q]  [-0.57721  +  Ln(l/q)  -  I  - ]  (13) 

n-1  n 

where  q  -  BA/( |rdotl*ri| 


The  median  range  of  acquisition  is  the  range  at  which  there  is  a 
50  percent  probability  that  visual  acquisition  has  occurred.  It  can  be 
written; 


r0*5 


1 

rdot|(Ln  2)/(6A)  +  1/ri 


(lA) 
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APPENDIX  B 


FILE  STRUCTURE  FOR  SEE  PROGRAMS 

The  SEE  programs  are  written  in  the  Pascal  computer  language  using  the 
Turbo  Pascal  compiler  from  Borland  International*  This  compiler  is  available 
for  a  variety  of  microcomputer  operating  systems,  including  CP/M  and  MSDOS.  A 
CP/M  version  is  described  below. 

Table  B.l  shows  the  files  that  exist  on  disk  for  the  CP/M  version  of 
SEEl.  All  text  files  are  indicated  by  ”Y'*  in  the  second  column  of  the  table. 
Text  files  can  be  modified  by  using  the  "non-document"  or  "data"  mode  of  any 
text  editor  (Caution:  Do  not  alter  these  files  using  the  word  processing  or 
"document"  mode  of  a  word  processing  program.  Such  a  mode  inserts  control 
codes  such  as  pag<>  breaks  and  soft  hyphens  into  the  data.  The  SEE  programs 
will  not  function  if  they  encounter  such  control  codes  in  a  data  file). 

TABLE  B.l 

FILES  FOUND  ON  THE  CP/M  DISK 


FILE  NAME 

TEXT 

DESCRIPTION 

SEEAC.B03 

Y 

Table  of  aircraft  types  (names  and  principal  areas). 

SEEMODS.B03 

Y 

Common  procedures  and  functions  for  SEEl  and  SEE2, 
including  EXPP,E2,GETVAR,  MODULEl,  M0DULE2,  M0DULE3, 
M0DULE4,  MO’jULES,  M0DULE6,  MODULE?,  READR  and  SGN. 

SEEVAR.B03 

Y 

Initialization  data  (variable  names,  nominal  values, 
and  conversion  factors  between  internal  and  external 
units). 

SEE1SET.B03 

Y 

Set-up  file  for  SEEl.  Contains  names  of  all  other  files 
used  by  SEEl.  Also  contains  printer  set-up  string* 

SEE2IN.B03 

Y 

Simulation  initialization  data  for  SEE2.  Contains 
definition  of  reference  and  subject  aircraft  pairs. 

SEE2SET.B03 

Y 

Set-up  file  for  SEE2.  Contains  names  of  all  other  files 
used  by  SEE2.  Also  contains  printer  set-up  string. 

SEE1.CMD 

N 

CP/M  command  file  for  SEEl  (compiled  version  of  SEEl). 

SEE1.B03 

Y 

SEEl  source  code  (Turbo  Pascal) 

SEE2.CMD 

N 

CP/M  command  file  for  SEE2  (compiled  version  of  SEE2). 

SEE2.B03 

Y 

SEE2  source  code  (Turbo  Pascal). 
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A  basic  set-up  file,  SEE1SET.B03,  Is  used  to  specify  file  names  for  all 
other  files  used  during  the  running  of  SEE1«  This  file  also  specifies  the 
printer  Initialization  string  to  be  sent  to  the  printer  prior  to  generation  of 
printed  output.  Figure  B-1  contains  the  listing  for  the  version  of 
SEEISET. B03  provided  on  the  start-up  disk.  The  corresponding  file  for  SEE2  is 
shown  in  Fig.  B-2. 

Aircraft  data  (names  and  principal  areas)  are  contained  In  file 
SEEAC.B03.  Figure  B-3  contains  a  listing  for  the  default  aircraft  file 
provided  on  the  start-up  disk. 

Disk  file  SEEVAR.B03  contains  the  names  of  the  Input /output  variables, 
their  nominal  values,  and  the  conversion  factors  used  to  convert  from  external 
to  Internal  units.  Internal  to  SEEl,  all  units  are  based  upon  nautical  miles, 
seconds,  and  radians.  Figure  B-4  provides  a  listing  -^f  the  default  file. 


SEEl  SET-UP  FILE  ;  Version  B03 
A:SEEVAR.B03  {VARIABLE  NAME  FILE} 

A:SEEAC.B03  {AIRCRAFT  TYPE  FILE) 

A: SEEl. OUT  {OUTPUT  DATA  FILE} 

27  91  50  119  -1  -1  {SIX  ASCII  CHAR  FOR  PRINTER  SET-UP:  -1  TO  SKIP} 


Figure  B-1  Input  data  file  SEEISET. B03 


SEE 2  SET-UP 

A;SEE2IN.B03 

A:SEEVAR.B03 

A:SEEAC.B03 

A;SEE2.0UT 

27  91  50 


FILE  Version  B03 

{SET-UP  DATA  FILE} 

{VARIABLE  NAME  FILE} 

{AIRCRAFT  TYPE  FILE} 

{OUTPUT  DATA  FILE} 

119  -1  -1  {SIX  ASCII  CHAR  FOR  PRINTER  SET-UP:  -1  TO  SKIP} 


Figure  B 


1 

PA28 

35. 

85, 

2 

B727 

400.  1900 

3 

B747 

1200.  5700 

4 

DC-9 

300.  1425 

5 

C421 

81. 

171 

6 

U-21 

127. 

267 

7 

F18 

50. 

280, 

2  Input  data  file  SEE2SET.B03 


260. 

3100. 

9300. 

2325. 

417. 

652. 

540. 


Figure  B-3  Aircraft  data  file  SEEAC.B03 
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NAME 

NCMINAL 

CONV 

EXT  UNITS 

BEI? 

17000.0 

1.0 

/STER-SEC 

BETl 

17000.0 

1.0 

/STER-SEC 

T1 

180.0 

1.0 

SEC 

T2 

12.0 

1.0 

SEC 

Dl 

0.0 

1.0 

NMI 

D2 

0.0 

1.0 

NMI 

R 

20.0 

1.0 

NMI 

DLIM 

1.0 

2.90883E-04 

ARC-MIN 

VSUB 

180.0 

2.7777E-04 

XT 

VSUB 

130.0 

2.7777E-04 

KT 

lACT 

1.0 

1.0 

AC  TYPE 

XANG 

180.0 

0.017453 

DEG 

PL 

1.0 

1.0 

LEFT  PILOT 

PR 

0.0 

1.0 

RIGHT  PILOT 

FOVL 

-120.0 

0.017453 

DEG 

FOVR 

90.0 

0.017453 

DEG 

ROOT 

300.0 

2.7777E-04 

KT 

AX 

35.0 

2.70872E-08 

SQ.  FT. 

AY 

85.0 

2.70872E-08 

SQ.  FT. 

AZ 

260.0 

2.70872E-08 

SQ.  FT. 

AREA 

35.0 

2.70872E-08 

SQ.  FT. 

THT 

0.0 

0.017453 

DEG 

THS 

0.0 

0.017453 

DEG 

RNGl 

3.78 

1.00 

NMI 

RNG2 

1.27 

1.00 

NMI 

BEFO 

17000. 

1.00 

/STER-SEC 

BEFl 

17000. 

1.00 

/STER-SEC 

PACQ 

0. 

1.00 

— 

Figure  B-4  Variable  initialization  file  SEEVAR.B03 
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APPENDIX  C 


SEE  LISTINGS 

This  appendix  provides  listings  for  the  SEEl  programs*  The  programs  were 
written  In  the  Pascal  computer  language  and  compiled  using  the  Turbo  Pascal 
compiler  from  Borland  International,  Incorporated*  Section  C* 1  contains  all 
procedures,  functions,  and  main  program  code  for  SEEl,  arranged  In 
alphabetical  order  of  name  (the  main  program  Is  listed  as  SEEl)*  Section  C*2 
contains  all  code  for  SEE2  with  the  exception  of  routines  that  have  been 
previously  listed  In  Section  C*l*  It  should  be  noted  that  In  Turbo  Pascal,  If 
a  procedure  B  calls  procedure  A,  then  A  must  be  listed  first  In  the  program 
listing*  Hence  the  actual  source  code  that  Is  compiled  will  not  have 
procedures  In  alphabetical  order*  The  source  code  also  uses  an  "include” 
directive  to  merge  the  common  procedures  In  file  MODULES  with  the  main  program 
listing  during  compilation* 

C* 1  Listing  for  SEEl 

PROCEDURE  DlsppCVar  pp,ppe,conv  :  ppty;  name  :  namety; 

Index, jndex  :  Integer); 

{DISPLAY  INPUT  VARIABLES  TO  CONSOLE} 

VAR  skip  ;  string [12];  l,j,li  :  integer; 

begin 

LowVldeo; 

Exunit8(pp,ppe,conv); 

For  1:“1  to  6  do 
begin 

For  j :«1  to  3  do 
begin 

ii:-6*(j-l)+i; 

If  (11<«16)  then 
begin 

WriteC  '.name  [11],*  •); 

If  (ppe[li]<10*0)  then  Write(ppe[il]; 10:2)  else  Write(ppe[il]; 10: 1); 
If  (ii*index)  or  (il»jndex)  then  WriteC*')  else  WriteC  ')» 
end; 
end; 

Writeln; 
end ; 

Writeln; 

WritelnC  *  ■  value  on  first  iteration'); 

NormVideo; 

end;  {of  PROCEDURE  Dispp} 

FUNCTION  Expp(x  :  real)  :  real; 

begin  {Exponentiation  with  overflow/underflow  protection) 

If  (x>85*0)  then  Expp:"l*0e+37 
else  If  (x<-85*0)  then  Expp:«l*0e-37 
else  Expp:«exp(x); 
end ; 
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FITNCTION  E2<x  :  real)  :  real; 

{Evaluate  exponential  integral  of  order  2  from  x  to  Infinity) 

Var 

df,fl,f2  :  real;  11, it  :  Integer; 

begin 

e2:-1.0E-37; 

If  (x<9.0)  then 
begin 

f 1 : “expp ( -X )+9. 57 72 l*x4x*Ln (x ) ; 
df x*x; 
f2:-df ; 
lt:-0; 

Repeat 

it:-lt+l;  li:-(lt+l)*(lt+l); 
df it*x*df/ii; 
f2:-f2+df ; 

Until  (ab8(df )<1.0E-09); 

e2:-fl+f2; 

end; 

end; 

PROCEDURE  ExunltsCVar  pp,ppe,conv  ;  ppty); 

{LOAD  PPE  WITH  VARIABLES  IN  EXTERNAL  UNITS) 

Var  1  ;  Integer; 
begin 

For  1:-1  to  28  do  ppe(l];«>pp[l]/conv[i] ; 
end;  {of  Exunits } 

FUNCTION  FpacqCVar  pp  :  ppty;  Var  tabl  :  tabty)  ;  real; 

Var  fov  :  real; 

begin  {COMPUTE  PROB  OF  VISUAL  FOR  CONDITIONS  IN  PP) 

Modulel(pp[ll],tabl,pp[18],pp[19],pp[20]); 
Module2(pp[9],pp[10],ppI12].pp[17],pp[22],pp[23]); 
Module3(pp[18],pp[19],pp[23],pp[21 j); 

Module4(pp[l],pp[2],pp[15],pp[16],pp[13],pp[14],pp[22],pp[26],pp[27],fov); 

Module5(pp[4],pp[6],pp[17],pp[25]); 

Module6(pp[3] ,pp[5],pp[17j,pp[24j); 

IF  (DEBU03)  THEN  BEGIN  WRITELN(Lst , 'INPUT  TO  MODULE  7:  AREA-’ ,pp  [21  ] :  10: 8, 
'  BEF0,BEFl-',pp[26]:9:0,pp[27]:9:C); 

WRITELN(Lst , '  DLIM-’ ,pp [8] : 10:8, ’  R-’,pp[7]:7:3, ’  RDOT-',pp[17]:9:5); 

WRITELNCLst,'  RNG1-’ ,pp [24] : 7: 3, ’  RNG2-’ ,pp [25] : 7: 3) ;  END; 

Module7(pp[21],pp[26],pp[27],pp[8],pp[7],pp[17],pp[24],pp[25],pp[28]); 

Fpacq :-pp[28] ; 

end;  {OF  FUNCTION  FPACQ) 

PROCEDURE  GetVar(Var  st  :  8trng4;  Var  index  :  Integer); 

{VALIDATE  VARIABLE  NAME  AND  GET  INDEX.  CR  RETURNS  INDEX— 1 } 

Var  1  :  integer; 

begin 

Repeat 
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Wrlte(’»');  Readln(8t);  index;-0; 

For  1:-1  to  Length(8t)  do  8t[i]:-UpCa8e(8tti]); 

While  (Length(8t)<4)  do  8t:“8t+’ 

If  (8t<>'  ')  then 

begin 
index 

For  i:-!  to  16  do  If  (name [1] -at)  then  index: “i; 

IF  (index<0)  then  Writeln(Chr(7), ’ILLEGAL  NAME  -  TRY  AGAIN’) 
end; 

Until  (index>"'0); 

end;  {  OF  GETVAR  MODIFIED  12/30/87} 

PROCEDURE  LatppCVar  pp,ppe,conv  :  ppty;  Var  name  :  namety); 

{PRINT  INPUT  VARIABLES  TO  PRINTER} 

Var  l,ll,j  :  Integer; 
begin 

Exunlt8(pp,ppe,conv) ; 

For  1:"1  to  6  do 
begin 

For  j:"l  to  3  do 
begin 

ii:-6*(j-l)+l; 

If  (il<-16)  then 
begin 

Write(OutF,'  ’,name[il],’  ’); 

If  (ppe [11 ]<10.0)  then  Wrlte(OutF,ppeIii] : 10:2)  elae  Wrlte(OutF»ppe[ii) : 10: 1) 
If  (ll*index)  or  (il»jndex)  then  Wrlte(OutF, **’)  else  Wrlte(OutF, ’  '); 
end; 
end; 

Writeln(OutF); 

end; 

Writeln(OutF); 

Wrlteln(OutF, '  *  ■  value  on  first  iteration'); 

Writeln(OutF); 

end;  {of  PROCEDURE  Lstpp} 

PROCEDURE  Module l(xlact  :  real;  tabl  :  tabty;  Var  ax,ay,az  :  real); 

Var  lact  :  Integer; 

{ASSIGN  TARGET  AIRCRAFT  AREAS  BY  INDEXING  TABl} 
begin 

lact : “Round (xiact ) ; 
ax: “tabl [lact , 1 ] *8qnml ; 
ay: “tabl [lact,2]*8qnml; 
az:“tabl[iact,3]*sqnmi; 

end;  {of  Procedure  Module  1  7/30/87} 

PROCEDURE  Module2(Var  vl,v2,xang,rdot,thl,th2  :  real); 

{COMPUTE  TH1,TH2  AND  RDOT  ASSUMNG  COLLISION  COURSE} 

Var 

dif ,cosmax,co8x,8inth,slnx,vco8max,vrat  :  real; 
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begin 

If  (Ab8(v2-vl)<0. 00001)  then  {USE  EQUAL  SPEED  FORMULAS} 

begin 

If  (xang“0. 0)  then  xangi“0« 00001;  {DISALLOW  DEGENERATE  CASE} 
thl : -xang/2-0. 5*pl*8gn(xang) ; 
th2:— thl; 

rdot : — 2*vl*ab8(8ln(xang/2 ) ) ; 
end 

else  {USE  UNEQUAL  SPEED  FORMULAS} 

begin 

vrat;»v2/vl;  co8x:-co8(xang); 

If  (xangOO.  0)  then 
begin 

th 1 : -arctan ( -vrat  *8in ( xang ) / ( 1 . 0-vrat*cosx ) ) ; 

If  (co8x>1.0/vrat)  then  thl:“thl-pl*sgn(thl) ; 
end 

elae  if  (vrat<1.0)  then  thl:"0.0 
elae  thl: -pi; 
t h2 : -p i +t h 1 -xang ; 

If  (ab8(th2)>pi)  then  th2:-th2-2*pl*8gn(th2); 

rdot ; ■“vl*co8(thl )-v2*co8 (th2) ; 

end; 

end;  {of  Module  2  Vera ion  3/11/87} 


PROCEDURE  Module3(ax,ay,tht  :  real;  Var  area  :  real); 

{COMPUTE  VISUAL  AREA  OF  TARGET.  USE  APPROXIMATE  CORRECTION  FOR 
HIDDEN  AREAS} 

Var 

axx,ayy  :  real; 
begin 

axx : ■ax*ab8 ( coa ( tht ) ) ; 
ayy :-ay*ab8(8in(tht) ) ; 

If  (axx>a3ry)  then  area:“axx+ayy/3  elae  area:“ayy+axx/3; 
end ; 


PROCEDURE  ModuleA(betO,betl,fovl,fovr,pl,pr,thl  :  real; 

Var  bef 0,bef l,fov  :  real); 

{DETERMINE  EFFECTIVE  VALUES  OF  BETA  BEFORE  AND  AFTER  TRANSITION} 

{FOV-0  IF  TARGET  NOT  IN  FIELD  OF  VIEW  OF  ANY  PILOT} 

begin 

bef0:“0.0;  befl:“0.0;  fov:“0.0; 

{IF  IN  LEFT  FOV,  ADD  CONTRIBUTION  OF  PILOT  IN  LEFT  SEAT} 

If  (thl>-fovl)  and  (thK-fovr)  then 
begin 

bef 0:-bef 0+pl*bet0; 
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bef 1 : “bef l+pl*bet 1 ; 

fov;«4ov+pl; 

end; 

{IF  IN  RIGHT  FOV,  ADD  CONTRIBUTION  OF  PILOT  IN  RIGHT  SEAT) 
If  (thl>"-fovr)  and  (thl<“-fovl)  then 
begin 

bef  0: “bef  0+pr*bet  0 ; 
bef 1 : "bef I4pr*bet 1 ; 
fov:“fov+pr; 
end; 

end;  (of  MbduleA  12/29/87} 


PROCEDURE  Module5(t2,d2,rdot  :  real;  Var  mg2  :  real); 

(USE  MODIFIED  TAU  TO  COMPUTE  RANGE  AT  WHICH  PACQ  IS  EVALUATED) 
begin 

rng2: -d2-rdot*t2; 
end; 


PROCEDURE  Module6(tl,dl ,rdot  :  real;  Var  rngl  :  real); 

(USE  MODIFIED  TAU  TO  COMPUTE  RANGE  AT  WHICH  BETA  TRANSITIONS) 
begin 

rngl:»dl-rdot*tl; 
end;  {OF  M0DULE6} 


PROCEDURE  Module7(area, bef 0, bef l,dlim,r,rdot, rngl, rng2  :  real;  Var  pacq  :  real); 
{COMPUTE  PACQ) 

Var 

absrd,pacq0,rllm,xlla,xpo8t ,xprlor,xx,xxl,xx2  :  real; 

begin 

abs  rd :  ■•abs  (  rdot  )+0. 0000 1 ; 
pacq0:*>0.0;  pacq:"0.0; 

rlim:-1.1284*8qrt(area)/dlim;  {RLIM  CANNOT  EXCEED  RESOLUTION  LIMIT  OF  EYE) 

If  (rliin>r)  then  rlim;-r;  {RLIM  CANNOT  EXCEED  VISUAL  RANGE} 

{ -  COMPUTE  INTEGRAL  OVER  THREE  POSSIBLE  TIME  INTERVALS) 

xlim:-e2(2. 996*rlim/r)/rlim; 
xxl:"e2(2.996*rngl/r)/rngl ; 
xx2:*e2(2.996*rng2/r )/rng2; 

If  (rlim>rngl)  then 

begin  {ADD  INTEGRAL  PRIOR  TO  TRANSITION) 

If  (rngl>rng2)  then  xx:“xx+bef 0*(xxl-Tcliin) 

elee  If  (rng2<rlim)  then  xx:”xx+bef 0*(xx2-xlim) ; 

end; 
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IF  (debug>4)  then 
begin 

WRITELN(L8t , ’  IN  MODULE?;  RLIM,RNG1,RNG2-' ,RLIM;8: 4,RNG1: 8: 4,RNG2; 8: 4) ; 
WRITELN(Lst, '  XLIM,m,XX2-',XLIM:8:4,XXl:8:4,XX2;8:4); 

WRITELN(Lst,'  XX  AT  BEPO-' ,XX; 10:4); 

end; 

If  (rngl>rng2)  then 

begin  {ADD  INTEGRAL  AFTER  TRANSITION} 

If  (rlim>rngl)  then  xx:“xx+bef l*(xx2-xxl) 

else  If  (rlim>rng2)  then  xx:*xx+bef l*(xx2-xlini); 

end; 

XX ; ■area*xx/ abs  rd ; 
pacq;*1.0-expp(-xx) ; 
end;  {OF  MODULE  7} 

PROCEDURE  0ut_con8t(Var  pp,ppe,conv  :  ppty;  Var  tabl  ;  tabty; 

Var  acname  :  acnamety;  Var  index, jndex  :  integer); 

{ -  PRINT  OUT  COMPUTED  VARIABLES  THAT  ARE  UNAFFECTED  BY  ITERATION  - 

Var  conll,conl7  :  Boolean; 
begin 

Modulel(pp[ll] ,tabl,pp[18],pp[19],pp[20]); 

Module2(pp[9],pptl0],pp[12],pp[17],pp[22],ppl23]); 

Module3(pp[18],pp[19i,pp[23],pp[21]); 

ExunitsCpp ,ppe ,conv) ; 

If  (index*il)  or  (Jndex*!!)  then  con!!:*FALSE  else  con!!;“TRUE; 
con !7: -TRUE; 

If  (index«9)  or  (Jndex*9)  or  (lndex«10)  or  (jndex"10)  then  con!7;»FALSE; 

If  (index«!2)  or  (jndex*!2)  then  con!7:*FALSE; 

If  con!!  then 
begin 

Writeln(OutF, '  |  Target  aircraft  is  acname [iact] ) ; 

Write(OutF,'  |  AX*' ,ppe[!8] :6: !, '  sq  ft  AY-' ,ppe[!9]:6: !, 

'  sq  ft  AZ-' ,ppeI20J :6: !, '  sq  ft'); 

If  con!7  then  Writeln(OutF, '  AREA*' ,ppe [2!]; 7: !, '  sq  ft') 

else  Writeln(OutF) ; 

end; 

If  con! 7  then 
begin 

Write(OutF,'  |  RDOT*' ,ppe(!7];7; 1, '  kt  THT-' ,ppe[22] ; 7: !, '  deg'); 

WriteCOutF,'  THS*',ppe[23j;7:!, '  deg'); 

end; 

Writeln(OutF) ; 

end;  {OF  PROCEDURE  0ut__con8t  } 

PROCEDURE  Printpp(Var  pp,ppe,conv  :  ppty;  name  :  namety); 

{PRINT-OUT  VARIABLES  IN  EXTERNAL  UNITS) 

Var  l,li,J  :  Integer; 
begin 

ExunitsCpp, ppe,conv) ; 
ii:*0; 
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For  1;*1  to  27  do 
begin 

If  (1<28)  then 
begin 

Wrlte(OutF,name[l] ,ppe [1] : 10: 2, '  ' ); 

11:-11+1; 

If  (11-5)  then  begin  Wrlteln(OutF);  ll:-0;  end; 
end; 

end; 

Wrlteln(OutF,name[28],ppe[28]: 10:4); 
end;  {of  PROCEDURE  Prlntpp} 

PROCEDURE  ReadR(Var  nltems  ;  Integer;  Var  buf  :  rbuff); 

LABEL  RETYPE; 

Var  l,len, valid  :  Integer;  ch  :  char;  s  :  strlng[255];  r  :  8trlngt32]; 
begin 

RETYPE:  nltems :-0;  r:-"; 

Readln(s);  s:-8+',';  len:-Length(8); 

For  1:-1  to  len  do 
begin 
ch:-s [1] ; 

If  (chO'  ')  and  (chO’,’)  then  r:-r+ch 
else 

begin 

If  (Length(r)>0)  then 
begin 

nltems : -nltems+1 ; 

Val(r,buf [nltems j , valid); 

If  (valldOO)  then 
begin 

Wrlteln(Chr(7) , '  Error  at  "',r,’"  (character  ’,valld:3, ’ )’ ) 
Wrlte( 'Retype  line'); 
goto  RETYPE; 

'-..j  < 

-  t 

t 

end; 

end; 

end; 

end;  {of  PROCEDURE  ReadR  Revised  l/lOa/86} 

Program  SEEl; 

{$R+} 

Const 

version- 'B03 ' ;  sqnml— 2. 70872E— 08;  nmlsec- 2.77777E— 04; 

Im-'  ';  nrow8max-128; 

Label  Ll.REDEF, START; 

Type 

namety  «  array [1.. 32]  of  string [4]; 
acnamety  -  array [ 1. . 24]  of  strlng[12]; 
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ppty  ■  array tl.. 32]  of  real; 
rbuff  “  array [1.. 32]  of  real; 

8trng4  -  strlngfA]; 

tabty  “  array[1..16,1..4]  of  real; 

Var 

InF,OutF,TempF  :  text; 
fovtX,xmln,xmax,xstep,xval  :  real; 

debug, i,lact,icol, index, ix,  j  ,joax,jiidex  :  Integer; 
long, menu, nacty,nrow8,nitein8  :  Integer; 
conv,nomval,pp,ppe  :  ppty;  buf  :  rbuff; 

Ibuff  :  array [1.. 32]  of  Integer; 
table  :  array [ 1. •nrowemax, 1* *8]  of  real; 

name  :  namety;  at  :  8trlng[A];  note, s, 88  :  8tring[255]; 
acname  :  acnamety;  ufov  :  Boolean; 

pavg,ppcol  :  array[1..8]  of  real; 
tabl  :  tabty; 

ACFNam,InFNam, OutFNam, printer, VFNam  :  strlng[l4]; 

{$!  SEEMODS.B03} 

PROCEDURE  0l8pp(Var  pp,ppe,conv  :  ppty;  name  :  namety; 
index, jndex  :  Integer); 

{DISPLAY  INPUT  VARIABLES  TO  CONSOLE) 

VAR  aklp  :  8trlng[12];  l,j,ll  :  integer; 

begin 

LowVldeo; 

Exunits (pp ,ppe , conv ) ; 

For  i:«l  to  6  do 
begin 

For  j!-l  to  3  do 
begin 

ii:-6*(j-l)+i; 

If  (li<-16)  then 
begin 

WriteC  ',name[li],'  '); 

If  (ppe [li ]<10.0)  then  Write(ppe{ii ] ; 10: 2)  else  Write(ppe[il]; 10: 1); 
If  (ll”lndex)  or  (ll*jndex)  then  WriteC*')  else  WriteC  '); 
end; 
end; 

Writeln; 

end; 

Writeln; 

WrltelnC  *  “  value  on  first  iteration'); 

NormVideo; 

end;  {of  PROCEDURE  Dlspp) 

PROCEDURE  Lstpp(Var  pp, ppe, conv  :  ppty;  Var  name  :  namety); 

{PRINT  INPUT  VARIABLES  TO  PRINTER) 

Var  i,ii,j  :  Integer; 
begin 

Exunits (pp , ppe , conv) ; 

For  i:*l  to  6  do 
begin 

For  j:“l  to  3  do 
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begin 

If  (11<-16)  then 
begin 

Write(OutF,'  ' ,name [ii] » '  '); 

If  (ppe[ii]<10.0)  then  Write(OutF,ppeIii]: 10:2)  else  Write(OutF,ppe [ii] : 10: 1) 
If  (ii-index)  or  (ii-jndex)  then  Write (OutF, •*’ )  else  Wrlte(OutF, '  '); 
end; 
end; 

Writeln(OutF); 

end; 

Writeln(OutF) ; 

Writeln(OutF, '  *  ■  value  on  first  iteration'); 

Writeln(OutF); 

end;  {of  PROCEDURE  Lstpp} 

PROCEDURE  Out_const(Var  pp,ppe,conv  :  ppty;  Var  tabl  :  tabty; 

Var  acname  :  acnamety;  Var  index, jndex  :  integer); 

{ -  PRINT  OUT  COMPUTED  VARIABLES  THAT  ARE  UNAFFECTED  BY  ITERATION  - } 

Var  con 11, coni 7  :  Boolean; 
begin 

Modulel(pp[ll],tabl,pp[18],pp[19],pp[20]); 

Module2(pp[9],pp[10],pp[12j,pp[17],pp[22J,pp[23]); 

Module3(pp[18],pp[19],pp[23],pp[21]); 

Exuni ts(pp,ppe,conv); 

If  (index-11)  or  (jndex-11)  then  coni 1: -FALSE  else  coni  1: -TRUE; 
coni  7; -TRUE; 

If  (index-9)  or  (jndex-9)  or  (index-lO)  or  (jndex-10)  then  conl7:-FALSE; 

If  (index-12)  or  (jndex-12)  then  conl7:-FALSE; 

If  conll  then 
begin 

Writeln(OutF, '  |  Target  aircraft  is  *, acname [ iact ]) ; 

Write(OutF, '  |  AX-' ,ppe[18]:6: 1, '  sq  ft  AY-' ,ppe [19] : 6: 1, 

'  sq  ft  AZ-' ,ppe [20] : 6: 1, '  sq  ft'); 

If  conl7  then  Writeln(OutF, '  AREA-',ppe[21]:7: 1, '  sq  ft') 

else  Wrlteln(OutF) ; 

end; 

If  con 17  then 
begin 

Write(OutF,'  \  RDOT-' ,ppe[17]:7; 1, '  kt  THT-’ ,ppe[22]:7: 1, '  deg'); 

Write(OutF,'  THS-' ,ppe[23] : 7; 1, '  deg’); 

end; 

Writeln(OutF); 

end; 

PROCEDURE  Printpp(Var  pp,ppe,conv  :  ppty;  name  :  namety); 

[PRINT-OUT  VARIABLES  IN  EXTERNAL  UNITS} 

Var  i,ii,j  :  integer; 
begin 

Exunlts(pp,ppe,conv) ; 
ll:-0; 

For  i:-l  to  27  do 
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begin 

If  (1<28)  then 
begin 

Write(OutF,naine  [i  ]  ,ppe[l] :  10:2,'  ' ) ; 

11 :“li+l ; 

If  (11-5)  then  begin  Writeln(OutF);  ii:*0;  end; 
end; 

end; 

Writeln(OutF,name[28] ,ppe [28] : 10:4) ; 


end;  {of  PROCEDURE  Printpp} 

begin  { - - - BEGIN  SEEl  MAIN  PROGRAM  -  -  } 


debug :-0; 

ClrScr;  GotoXY(15,l); 

WrlteCSEEl  VISUAL  ACQUISITION  ANALYSIS  -  Version  '.version); 
GotoXY(30,5);  Wrlte( 'written  for  ');  GotoXY(19,7) ; 

UrlteC'The  Federal  Aviation  Administration'); 

GotoXY(30,8);  Write (' (AT-240) ') ;  GotoXY(34, 10);  Write('by'); 
GotoXY(23,12);  Wrlte( 'M.I.T.  Uncoln  Laboratory');  GotoXY(17,14); 

WrlteC '244  Wood  St.,  Lexington,  MA  02173-0073');  GotoXY(l,23); 

WriteC 'CR  to  continue»');  Readln(st); 

If  (8t“'menu')  then 
begin 

Wrlte( 'Type  debug»');  Read  In  (  debug  ) ; 
end; 

InFNam:-'SEElSET.B03';  WrltelnC  1  INIT  FILE  NAME:  '.InFNam); 
Assign(InF,InFNam) ;  Reset(lnF); 

ReadlndnF.s);  Writeln(  ’  (  INIT  FILE  TITLE.-'.s); 

Readln(InF, VFNam) ;  Readln(InF,ACFNam) ;  Readln(InF,0utFNam) ; 

{  -  READ  PRINTER  INITIALIZATION  STRING  AND  INITIALIZE  PRINTER  } 

Readln(InF,lbuff [ 1 ] ,lbuff [2] ,lbuf f [3] ,lbuff [4] ,lbuff [5] ,lbuff [6] ) ; 
printer:-'';  For  i:-l  to  6  do 

If  (lbuffflj>0)  then  printer ;-prlnter+Chr(lbuff [1 J ) ; 

If  (printerO'')  then  WrlteCLst, printer); 

{  -  READ  TABLE  OF  AIRCRAFT  CHARACTERISTICS  -  } 

Asslgn(TempF, ACFNam) ;  Reset(TempF); 
nacty:-0; 

While  not  Eof(TempF)  do 
begin 

nacty:-nacty+l; 

Readln(TempF, i .acname [nacty ] , tabi [nacty, 1  ] , tabl [nacty , 2 ] , tabl [nacty , 3 ] ) ; 
end; 

Close(TempF); 

{  -  read  in  variable  NAMES  AND  NOMINAL  VALUES  -  } 

AssignCTempF, VFNam) ;  Reset(TempF) ; 

Readln(TempF);  (Skip  label  line} 

For  i:-l  to  28  do 

Readln(TempF,name [1 1 .nomval [l] ,conv[i] ) ; 

Close (TempF ) ; 

WriteC 'Type  1  for  long  output  (CR  for  summary  table  only)»'); 
long:-0;  Readln(long) ; 
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START: 

For  1:«1  to  28  do 
begin 

ppe[l  ]  :«noinval  (1 1 ; 

pp[l):“nomvalIi]*conv[i] ;  {Initialize  variables  to  nominal  values} 
end ; 

{  -  SET  INPUT  CONDITIONS  -  } 

REDEP  : 

Repeat 

Write(’Type  name  of  row  variable  ');  GetVar(st, index); 

Until  (lndex>0}; 

WriteC’Type  xmln,xmax,xstep  for  *  ,name[lndex] ,  *»' ); 

ReadR(nltems,buf ) ;  xmln : -buf [ 1 ] ;  xmax : -buf [ 2 ] ;  xstepr-buf [3] ; 
nrows : “Round^ (xmax-xmln)/xstep+0. 5) ; 

If  (nrows >nrowsmax)  or  (nrows<“0)  then 
begin 

Writeln(Chr(7), ’»»  NO.  OF  ROWS  MUST  BE  BETWEEN  1  AND  ’ , nrowsmax : 5 ) ; 
Writeln('»»  PLEASE  INPUT  ROW  DATA  AGAIN.*); 

Goto  REDEF; 
end; 

pp [index] :“xnln*convIlndex] ; 

Repeat 

WrlteC'Type  name  of  column  variable  ');  GetVar(st,Jndex) ; 

Until  (jndex>0); 

Writeln( 'Type  up  to  8  values  for  *,  namel jndex] , '  »’); 

Readr(jmax,buf ) ;  For  jS“l  to  jmax  do  ppcoltj ]:“buf [j J; 
pp [Jndex]: "ppcol [ 1 ] *conv [jndex ] ; 

-  SELECT  TARGET  AIRCRAFT  TYPE  -  } 

WritelnC 'TABLE  OF  DEFINED  TARGET  AIRCRAFT  TYPES  :’); 

WrltelnC  lACT  NAME  AX  AY  AZ'); 

For  1:«1  to  nacty  do 

WRITELNC  ',i:3,'  ' ,acname[l ] , ’  ',tabl[i,l]:6:0,tabl[i,2]:8:0, 

tabl[l,3]:8:0); 

Wrlte('Type  lACT  for  target  aircraft  (CR  to  define  new  type)»'); 
lact:“-l;  Readln(lact) ;  pp[ll]:“1.0*iact; 

If  (iact<“0)  then 
begin 

nacty :*nacty+l;  lact: -nacty;  pp[ 11 ] ;-1.0*nacty ; 

WriteCType  symbolic  name  for  type  (12  char  max»'); 

Readln(acname[lact] ); 

While  (Length(acname [lact ] )<12)  do  acname[lact]:-acname[iact]+'  '; 

Writeln( 'Type  AX.AY.AZ  (in  sq  ft)»’); 

ReadR(nltems , buf ) ; 

tabl [lact, 1 ] : -buf [ 1 ] ;  tabl [lact ,2]:-buf [2] ;  tab 1 [lact ,3] : -buf [3] ; 
end; 

LI  : 

{  -  DISPLAY  SET-UP  TO  CONSOLE  FOR  APPROVAL  -  } 

ClrScr;  GotoXY( 10, 1) ;  Write ( '*-*-*-*-*-*-*  SEEl  ANALYSIS  -  Version  version); 
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Wrlteln( '  *-*-*-*-*-*-* • ) ; 

WritelnCROW  VARIABLE  ;  ');  LowVideo; 

WrltelnC  name [index] , '  to  be  varied  from  ',xmln:8;l,'  to* ,xmax:9: 1, 

'  in  steps  of  ',xstep:8:2);  NormVldeo; 

WrltelnC 'COLUMN  VARIABLE  2  ');  LowVideo; 

WriteC  ' ,name[jndex] , '  values  :*); 

For  1:"1  to  Jmax  do  Write(ppcol[l]:8:2);  NormVldeo;  Writeln; 

WrltelnC 'TARGET  AC  :  WITH  AREAS  CAX.AY.AZ)  :  ');  LowVideo; 

WrlteC '  ',acname[lact] , '  '); 

WrltelnCtabl [iact , 1 ] : 7: l,tabl [lact ,2] :7: l,tabi [lact,3] : 7:1);  NormVldeo; 

Writeln; 

pp t Index ] :“xmin*conv[ index] ;  pp[ jndex]:-ppcol [ 1 ]*conv[jndex] ; 

WrltelnC 'CURRENT  VALUES  OF  INPUT  VARIABLES;’);  Writeln; 

DlsppCpPfPpe, conv, name, index, jndex);  {Display  input  variables} 

WrlteC 'Type  name  of  any  variable  to  be  changed  CCR  if  none)'); 

GetVarCst,!); 

If  Ci>0)  then 
begin 

WrltelnC 'Currently  ' ,name [i ] , *■' ,ppe[i] : 10:2); 

WriteC  *Type  new  value»');  ReadRCnitems,buf ) ;  x:“buf[l]; 
ppe [ i ] ; -X ;  pp I 1 ] : ■x*conv [ 1 ] ; 

If  Cl"!!)  then  iact :->RoundCpp [!!]);  goto  L!;  end; 

[ - SET-UP  COMPLETE.  RUN  ANALYSIS  -  -  } 

WrltelnC 'MENU  ;');  WrltelnC  CR  -  RUN  ANALYSIS’); 

WrltelnC  !  -  REDEFINE  INPUT  CONDITIONS’); 

WrlteC»');  12—44;  ReadlnCi); 

If  ClO-44)  then  goto  REDEF; 

WrlteC'Type  note  CCR  if  none)»');  ReadlnCnote) ; 

{ - WRITE  HEADING  INFO  TO  OUTPUT  FILE  OUTF - } 

AsslgnCOutF,OutFNam) ;  RewrlteCOutF) ; 

WritelnCOutF.lm, 'SEE!  VISUAL  ACQUISITION  ANALYSIS'); 

If  CbengthCnote)>0)  then  WrltelnCOutF,lm, '  |  NOTE  :',note); 

WrltelnCOutF,lm, '  SEE!. PAS  CVersion  '.version,')'); 

WritelnCOutF.lm,'  InFNam-  ' ,InFNam; 17, '  VFNam-  ',VFNam); 

WritelnCOutF.lm,'  ACFNam*  '  ,ACFNam;  !7. '  OutFNam*  '-.OutFNam); 

iact: >RoundCpp [!!]); 

WrltelnCOutF) ; 

WritelnCOutF.lm, 'INITIAL  VALUES  OF  INPUT  VARIABLES:'); 

Ls tppCpp, ppe, conv, name ) ; 

WritelnCOutF.lm, 'ROW  VARIABLE:'); 

WrltelnCOutF, Im, Im, name[lndexl, '  is  varied  f rom' .xrnln: 9: 2, '  to', 
xmax:9:!,'  in  steps  of ' ,xstep:8:2); 

WritelnCOutF.lm, 'COLUMN  VARIABLE  :’); 

WriteCOutF,lm,lm,name[jndex], '  assumes  values  '); 

For  1:"!  to  jmax  do 
begin 

If  Cl-5)  then 

begin  WrltelnCOutF);  WrlteCOutF.lm.lm.lm.lm.lm, '  ');  end; 
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Write (OutF.ppcol [1 ] : 8; 2) ; 
end; 

Wrlteln(OutF) ; 

Out  con8t(pp,ppe,conv,tabl,acname, index, jndex) ; 

{  - - START  COLUMN  ITERATION - } 

For  lcol:"l  to  jmax  do 
begin 

WrltelnC  NOW  CtWPUTING  WITH  ' ,name(jndex] , .ppcol [icol] : 8: 2) ; 
pp  t J  ndex ] : “ppcol [ Icol ] *conv [ jndex ] ;  pavg ( icol ] : -0. 0 ; 

If  (long“l)  then 
begin 

Wrlteln(OutF); 

Writeln(OutF,lm,'SEEl  ANALYSIS  :  PACQ  VS.  ' ,name[lndex] , ' , 
name [Jndex] , , ppcol [Icol ] :8:2); 

Writeln(OutF,lm, '  name [index] , '  RDOT  AREA  THT  THS', 

’  RNGl  BEFl  PACQ*); 
end; 

{ - START  ROW  ITERATION - } 

lx;“0;  xval;“xmin-xatep; 

While  (xval<xmax)  do 

begin 

lx:“lx+l ; 

xval :"xval+x8tep;  ppe[ index ]:“xval;  pp [index] :“xval*conv[ index ] ; 

Modulel(pp( 11 ],tabl,pp[18],pp[19],pp[20]); 

Module2(pp[9],pp[10],pp[ 12],pp[17],ppl22],pp[23]); 

Module3(pp(18],pp[19],pp[23] ,pp[21]); 

Module4(pp[l],pp[2],pp(15],pp[16],pp[13],pp[ J4],pp[22],pp[26],ppI27],fov); 
Module5(pp(4] ,pp(6],pp[17l,pp[25] ); 

Module6(pp [ 3 ] , pp [ 5  ]  ,pp [ 1 7 ] , pp [ 24 ] ) ; 

IF  (DEBUO-2)  THEN  BEGIN  WRITELN(OutF,  ’  INPUT  TO  MODULE  7:  AREA-’ ,pp [21  ] ;  9:  6, 
’  BEF0-',pp(26];7:2, ’  BEF1-* ,pp [27 J: 7; 2) ; 

WRITELN(OutF, *  DLIM-' ,pp[8] : 7:2, •  R-’ ,pp[ 7 ] : 7: 2, ’  RDOT-' ,pp[ 17] : 7; 3) ;  END 
Module7(pp[21],pp[26],pp[27],pp[8],pp(7],pp[17],pp[24],pp[25],pp[28]); 
pavg [icol ] : -pavg [ icol ]+pp [28] ; 

Exunit8(pp,ppe,conv) ; 

If  (long-1)  then  Wrxteln(OutF,lm,ppe[indexl: 9; l,ppe[ 17 ] : 9; l,ppe [21 ] : 9: 1, 
ppe [22] : 7: l,ppe[23] : 7: l,ppe(24] : 7:2,ppe[27] : 9;0,ppe [28 j : 9: 4); 

If  (lx<32)  and  (fov>0.0)  then  tableIix,icol];-ppe[28] 
else  if  (fov-0.0)  then  table[ix, icol]:— 1.0; 
end;  {of  index  block] 

{ - END  ROW  ITERATION - } 

pavg [icol ] : -pavg [icol ] /nrows ; 
end;  {of  icol  block] 

(  -  end  COLUMN  ITERATION  - ] 

{ - - - OUTPUT  PACQ  TABLE - } 

Wrlteln(OutF, 'PACQ  TABLE  :  '); 

Writeln(OutF); 

Writeln(OutF,lm,'  TABLE  OF  PACQ  VALUES’); 

Wrlteln(OutF) ;  Writeln(OutF,lm, '  ', name [jndex] ) ; 


Wrlte(OutF,lni,  '  ');  ^  ^ 

For  1:-1  to  jmax  do  Write(OutF,ppcol [1 ] : 9: 2) ;  Writeln(OutF) ; 

Write(OutF,lin,  ’  '  .name  [Index] ,  '  '); 

For  1:“I  to  jmax  do  Wrlte(OutF,'  -  *); 

Wrlteln(OutF); 
ufov: "FALSE; 

For  1:“1  to  ix  do 
begin 

x:-xmin+(i-l)*x8tep;  Write(OutF,lo,x: 8; 1); 

For  j:“l  to  jmax  do 

If  (table [l.j ]>-0. 0005)  then  Wrlte(OutF, table[l, j ] : 9; 4) 
else  if  (table [1 ,j ]"-1.0)  then 

begin  Wrlte(OutF,'  xfov  *);  ufov:“TRUE;  end 
else  Wrlte(OutF,'  *  *); 

Wrlteln(OutF); 

end; 

pp [ Index ] : "nomval [ Index ] *conv [ Index ] ; 
pp[ jndex] :"nomval Ijndex]*conv [ jndex] ; 

Wrlteln(OutF) ;  Wrlte(OutF,lm, '  avg.'); 

For  1:-1  to  jmax  do  Wrlte(OutF,pavg[l]:9:4);  Wrlteln(OutF); 

If  ufov  then  begin  Writeln(OutF) ;  Wrlteln(OutF,lm,lm,lm,lm, 

'(  Note  :  xfov  ■  target  not  within  fleld-of-vlew) ' );  end; 

Wrlteln(OutF) ; 

Wrlteln(OutF,lm, '*-*-*-*  END  OF  ANALYSIS  *-*-*-*•); 

Close(OutF) ; 

Writeln( 'OUTPUT  STORED  IN  DISK  FILE  [ ' .OutFNam, ’ ] ' ) ; 

{ - SELECT  MODE  FOR  OUTPUT - } 

Wrlteln( 'MENU: ' ) ;  WritelnC  1  -  print  results'); 

Wrlteln('  2  ■  display  results');  Wrlteln( '  CR  ■  continue'); 
Wrlte('»');  menu:"0;  Readln(menu) ; 

If  (menu>0)  then 
begin 
ClrScr ; 

If  (menu"l)  then  AuxOutPtr:"LstOutPtr  else  AuxOutPtr:"ConOutPtr ; 
A8slgn(0utF, OutFNam) ;  Re8et(0utF); 

While  not  Eof(OutF)  do 
begin 

Readln(0utF,8 ) ; 

{If  (Length(s)>85)  then  WRITELNC  LENGTH(S)-' ,Length(8):4); } 

ss: "Copy(8 , 1, 10) ; 

If  (ss-'PACQ  TABLE')  and  (menu-2)  then 

begin  WrlteCCR  to  view  PACQ  table»');  Readln;  ClrScr;  end 
else 

begin 

If  (menu-2)  and  (Copy(8, 1,5)-'  ')  then  s: -Copy (8,6, 250); 

Wrlteln(Aux,8 ) ; 
end ; 

end ; 

Clo8e(0utF) ; 
end ; 
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-SELECT  TERMINATION  OPTIONS  -  } 


Write( 'CR  to  contlnue»');  Readln; 

{ - 

WrltelnC’MENU 

WrltelnC  1  -  NEW  ANALYSIS  (KEEP  PREVIOUS  MODS)’)*, 

WritelnC  2  -  NEW  ANALYSIS  (REINITIALIZE  ALL  INPUTS)’); 

WrltelnC  CR  -  QUIT’); 

Wrlte(’»’);  me.iu:-0;  Readln(menu) ; 

If  (inenu“l)  then  goto  LI  else  if  (menu“2)  then  goto  START; 
Clo8e(0utF) ; 

WritelnC  OUTPUT  SAVED  IN  ’.OutFNam); 
end. 


C. 2  Listing  for  SEE2 

This  section  provides  a  listing  of  the  Pascal  code  for  SEE2.  Each 
procedure  is  shown  seprately  in  alphabetical  order  of  name.  Those  procedures 
that  are  identical  to  those  used  in  SEEl  are  not  included  (they  are  provided 
in  section  C.2). 


PROCEDURE  Dispp2(Var  pp,pp3,pp4,conv:  ppty;  name:  namety; 
xndex , yndex :  Intege r ) ; 

{DISPLAY  INPUT  VARIABLES  TO  CONSOLE} 

VAR  stx:  string [255]; 
begin 

WrltelnC  -  REFERENCE  PAIR  -  DUMMY  SUBJECT’); 

WritelnC  ACl  AC2  AC3  AC4’); 

For  i:*l  to  16  do 
begin 

8tx:*name  [i  ]  +  ’  -  -  -  - 

If  (i“xndex)  then  Writeln(8tx, '  ROW  VARIABLE’) 

else  if  (i*yndex)  then  Writeln(8tx, ’  FLOATING  VARIABLE’) 

else  Writeln(StDl8p(i,pp,pp3,pp4,conv,name)) ; 

end; 

Wrlteln; 

end;  {of  PROCEDURE  Dlspp2  Revised  7/30a/87} 


PROCEDURE  Make_Conslstent(Var  pp,pp3,pp4;  ppty); 

{ -  MAKE  V1,V2,XANG,AND  R  CONSISTENT  BETWEEN  PAIRS) 

begin 

pp3[9]:-pp4[10];  pp3[10]:-pp4I9] ;  pp3[7 ] ; »pp4[7] ;  pp[7] :-pp4[7 ] ; 

pp3[12]:-pp4[12]; 

end; 
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FUNCTION  StDisp (Index:  Integer;  Var  pp,pp3,pp4,conv:  ppty; 

Var  name:  namety):  8tr255; 

Var  vout:  array[1..4]  of  real;  1, format:  Integer;  8t,8x:  8trlng[255]; 

{  RETURN  PP  VARIABLE  PRINT-OUT  LINE  FOR  4  AC  } 

begin 

voutt 1 ]s*pp[ index ]/conv[ index] ; 
vout [2] :-pp [index] /conv [index] ; 
vout [3] :“pp3[ index ]/conv[index] ; 
vout [4  J : “pp4 [index ] /conv [index ] ; 
st:-name[index]+'  '; 
format:-! ; 

Caae  index  of 

1.. 2  :  format :-0;  5*. 6  :  format: -2; 

8  :  format:“2;  11  :  format:-©; 

15. . 20  :  format:-0; 
end ; 

For  i:-l  to  4  do 
begin 

If  (format-1)  then  Str(vout [1] : 10: l,8x) 
elae  If  (format-0)  then  Str(vout [1]: 10:0,8x) 
el8e  Str(vout[l]: 10:2,8x); 

8t :-8t+8x; 
end; 

StDisp:-8t; 

end;  (of  Function  StLisp} 


PROCEDURE  Outpp(Var  pp,pp3,pp4,ppe,conv:  ppty;  Var  name:  namety; 
xndex,yndex:  Integer); 

{Print  input  variables  to  OutF} 

Var  8tx  :  8tring[255]; 
begin 

Wrlteln(OutF, '  -  REFERENCE  PAIR  -  DUMMY  SUBJECT’); 

Writeln(OutF, ’  ACl  AC2  AC3  AC4'); 

For  1;-1  to  16  do 
begin 

8tx:-name[l]  +  '  -  -  -  -  »; 

If  (1-xndex)  then  Writeln(0utF,8tx, ’  ROW  VARIABLE') 

elae  if  (1-yndex)  then  Wrlteln(0utF,8tx, ’  FLOATING  VARIABLE') 

else  Wrlteln(OutF, StDi8p(i ,pp ,pp3,pp4.conv.name ) ) ; 

end; 

Writeln(OutF) ; 

end;  {of  PROCEDURE  Outpp} 

Program  SEE2;  {Runs  from  default  disk) 

Const 

version- 'B03 ‘ ;  sqnmi-2. 70872E-08;  nmi8ec-2. 77777E-04; 
default-0.44444;  debug-0; 
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Label  ENDITT,L1,L2,REDEF, START; 

Type 

str255  ■  8trlng[255]: 

namety  -  array [1.. 32]  of  string [4]; 

ppty  ■  array [1*. 32]  of  real; 

rbuff  ■  array [1.. 32]  of  real; 

strngA  -  string [4];  strng255  -  string [255]; 

tabty  -  array [ 1. « 16, 1. . 4]  of  real; 

Var 

TFlle,OutF,InF  :  text; 

area,ax,ay,az,befO,bef I,bet0,betl,dfdy,dlim,dy,dl,d2  :  real; 
fnew,fold,fovl,fovr,f8top,pacq,pacqref ,pl,pr  :  real; 
r,rdot,rngl,rng2,dy8top,ti,t2,thl,th2,vl,v2  :  real; 
x,xang,xniln,xmax,xsmax,x8tep,xval,y,ye,yflr8t,yll  :  real; 
ynew,yold,y8tep,yul ,y0  :  real; 
i,ll,lac,lact,index,lLST,lx,j ,jmax  :  integer; 
nacty,nitems,nitt,n8top,xndex,yndex  :  Integer; 
conv,nomval,pp,pp3,pp4,ppe  :  ppty;  buf  :  rbuff; 

table  :  array [1* * 32, 1. .8]  of  real; 
name  :  namety;  st  :  8tring[4]; 

ACFNam,InFNam,OutFNam,SetFNam,VFNam  :  string [12]; 

acname  :  array[l.«16]  of  string[12];  note, s, title  :  8trng255; 

acid  :  array[1..16]  of  Integer; 

pacqac,pcol  :  array[l*«8]  of  real; 

tabl  :  tabty; 

{$!  SEEM0DS.B03} 

PROCEDURE  Make  Consistent (Var  pp>pp3>pp4  ;  ppty); 

{ - MAKE  V1,V2,XANG,AND  R  CONSISTENT  BETWEEN  PAIRS} 

begin 

pp3[9]:-pp4[10];  pp3[10]:-pp4[9] ;  pp3[7]:-pp4[7];  pp[7]:-pp4[7]; 

pp3[12] :-pp4[12] ; 

end; 

FUNCTION  StDl8p(lndex  :  Integer;  Var  pp,pp3,pp4,conv  ;  ppty; 

Var  name  :  namety)  :  8tr255; 

Var  vout  :  array[l.<4]  of  real;  1, format  ;  Integer;  st,8x  :  string[255]; 

{  RETURN  PP  VARIABLE  PRINT-OUT  LINE  FOR  4  AC  > 

begin 

vout [ 1 ] : -pp [ index ] / conv [ index ] ; 
vout [ 2 ] : -pp [ index ] /conv [ index ] ; 
vout [ 3 ] : “PP  3 [ index ] / conv [ index ] ; 
vout [4 ] :»pp 4 [ index ]/conv[ index] ; 

St: “name [index ]+'  '; 

format :“1; 

Case  index  of 

1.. 2  :  format :“0;  5.. 6  ;  format: “2; 

8  :  format:“2;  11  :  format:“0; 

15..  20  :  format :“0; 
end; 
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For  1:*1  to  4  do 
begin 

If  (fonnat-1)  then  Str(vout [i] : 10: l.sx) 

else  If  (format-O)  then  Str(vout [1] : 10:0,sx) 

else  Str(vout[l]: 10:2,sx): 

st:“st+sx; 

end; 

St’?t8p:“st; 

e  iC ,  {  of  Function  StDlsp } 

PROCEDURE  Dlspp2(Var  pp,pp3,pp4,conv  :  ppty;  name  :  namety; 
xndex.yndex  :  Integer); 

{DISPLAY  SEE2  INPUT  VARIABLES  TO  CONSOLE  } 

VAR  stx  :  string [255]; 
begin 

WrltelnC  -  REFERENCE  PAIR  -  DUMMY  SUBJECT'); 

WrltelnC  ACl  AC2  AC3  AC4'); 

For  1:-1  to  16  do 
begin 

stx:  “name  [1]+'  -  -  -  - '; 

If  (l“xndex)  then  WrltelnCstx, '  ROW  VARIABLE') 

else  If  (l*yndex)  then  WrltelnCstx,'  FLOATING  VARIABLE') 

else  Wrlteln(StDl8p(l,pp,pp3,pp4,conv,name)) ; 

end ; 

Wrlteln; 

end;  {of  PROCEDURE  Dlspp2  Revised  7/30a/87} 

PROCEDURE  OutppCVar  pp,pp3,pp4,ppe,conv  :  ppty;  Var  name  :  namety; 
xndex,yndex:  integer); 

{PRINT  INPUT  VARIABLES  TO  OUTPUT  FILE,  OUTP) 

Var  stx  :  8trlng[255]; 
begin 

WrltelnCOutF,'  -  REFERENCE  PAIR  -  DUMMY  SUBJECT'); 

WrltelnCOutF, '  ACl  AC2  AC3  AC4'); 

For  !:■!  to  16  do 
begin 

8tx:“n£me[l]+'  -  -  -  - '; 

If  (l*xndex)  then  WrltelnCOutF, stx, '  ROW  VARIABLE') 

else  If  Cl"yndex)  then  WrltelnCOutF, stx, '  FLOATING  VARIABLE') 

else  WrltelnCOutF, StDl8pCl,pp,pp3,pp4,conv, name)); 

end; 

WrltelnCOutF); 

end;  {of  PROCEDURE  Outpp} 

begin  { - - - SEE2  MAIN  PROGRAM  -  -  } 

WrltelnC'*-*-*-*  SEE2  VISUAL  ACQUISITION  ANALYSIS  CVerslon  ',ver8lon,')  *-*-*-*') 
{ - READ  DATA  FROM  INITIALIZATION  FILE - } 

InPNam:-'SEE2SET.B03';  WrltelnC  |  INIT  FILE  HEADING  :',InFNam); 
A8slgnCInF,InFNam);  ResetCinF); 

ReadlnCInF,8);  WrltelnC  |  INIT  FILE  TITLE  :  [',8,']'); 

ReadlnCInF.SetFNam);  ReadlnCInF.VFNam);  ReadlnCInF,ACFNam) ; 

ReadlnCInF,OutFNam) ;  WrlteC  I  SetFNam*' ,SetFNam) ; 

WrltelnC  VFNam  -',VFNam);  WrlteC  |  ACFNam-' ,ACFNam); 

WrltelnC  OutFNam-' ,0utFNam) ; 
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Readln(InF,acid[l],acid[2],acid[3],acld[4],acldI5I,acld[6]); 
a:-";  For  1:-1  to  6  do  If  (acid[l]>0)  then  8:-8+Chr(acidIi] ) ; 

Wrlte(L8t,8): 

A88lgn(0utF,(}utFNam) ;  Rewrite (OutF) ; 

{ - read  in  variable  names  and  nominal  values - } 

A88lgn(TFile,VFNam) ;  ResetCTFlle); 

Readln(TFlle):  {Skip  label  line} 

For  1:“1  to  28  do 

Readln(TFlle,name [1 ] ,nomval[l] ,conv[l] ) ; 

Clo8e(TFlle); 

{ - - read  in  data  FROM  SET-UP  FILE - } 

Wrlteln( 'Type  name  of  file  containing  Initialization  data.'); 

Wrlte('(CR  for  '.SetFNam,')  »'); 

ReadlnCa);  If  (Length(8)>0)  then  SetFNamt^e; 

A88lgn(TFlle, SetFNam) ;  ResetCTFlle ) ; 
yflrst:  "default;  3rndex:"9; 

Readln(TFlle,s);  WrltelnCs); 

ReadlnCTFlle, title) ;  Wrlteln(tltle) ; 

Readln(TFlle,st ,xndex,xmln,xmax,x8tep) ; 

ReadlnCTFlle ,yll ,yul .dystop ) ; 

ReadlnCTFlle); 

For  1:"1  to  16  do 
begin 

ReadlnCTFlle, 8t, X, pp[ 1], pp3 [1], pp4[l ] ); 
end; 

For  1:«1  to  16  do 
begin 

pp[ll !“Pplll*conv[l] ;  pp3[l] :“pp3[l]*conv(l] ;  pp4tll :“pp4[l]*conv[l ] ; 
end; 

CloseCTFlle); 

{ -  READ  TABLE  OF  AIRCRAFT  CHARACTERISTICS  -  } 

AsslgnCTFlle, ACFNam) ;  ResetCTFlle) ; 
nacty:-0; 

While  not  EofCTFlle)  do 
begin 

nacty : “nacty+l ; 

ReadlnCTFlle , acid [nacty ] , acname [nacty ] , tabl [nacty , 1 ] , tab 1 [nacty , 2 ] , 

tabl [nacty, 3] ) ; 

end; 

CloseCTFlle); 

Make__Consl8tentCpp,Pp3,pp4); 

{ - PRINT-OUT  SET-UP  DATA - } 

WrltelnC 'SET-UP  DATA  FROM  FILE  ', SetFNam, ':') ; 
Dlspp2Cpp,pp3,pp4,conv,name,xndex,yndex) ; 

WrltelnC 'FLOATING  VARIABLE:  ',name[yndex], '  for  AC4  will  vary  to  force  equivalency'); 
WrltelnC 'ROW  VARIABLE:  ',name[xndex], '  stepped  f rom' ,xmln:8: 1, '  to' ,xmax: 9: 1 ) ; 
WrltelnC  In  steps  of '  ,x8tep:9:2); 

REDEF  : 

{ -  ALTER  EXISTING  PP  VALUES  -  } 
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LI  : 

Write ('Type  name  of  any  variable  to  be  changed  (CR  if  none)'); 

GetVarCst.l); 

If  (i“xndex)  or  (i“yndex)  then 

Writeln<Chr(7), '  WARNING  :  THIS  IS  A  PROGRAM-CONTROLLED  VARIABLE!'); 

If  (i>0)  then 
begin 

WrlteCType  AC  to  change  (l-4):»');  Readln(lac); 

Write  ('Type  new  value»');  ReadR(nltems,buf ) ;  x:“buf[l]; 

If  (iac<3)  then  pp[i]:“x*conv[i] 

else  If  (lac=3)  then  pp3[il :-x*convIi] 
else  If  (lac“4)  then  pp4[l] :-x*conv[l ] ; 

Make__Conslstent  (pp ,  pp3,  pp4) ; 

Wrlteln( 'REVISED  SET-UP  VARIABLES  :'); 
Dlspp2(pp,pp3,pp4,conv,name,xndex>yndex) ; 
goto  LI; 
end; 

{ - WRITE  HEADER  INFO  TO  OUTPUT  FILE,  OUTF - } 

Wrlte( 'Type  note  (CR  if  none)»');  Readln(note) ; 

Wrlteln(OutF, '*-*-*-*  SEE2  VISUAL  ACQUISITION  ANALYSIS  *-*-*-**); 

If  (Length(note)>0)  then  Writeln(OutF, 'NOTE: ' .note) ; 

Wrlteln(OutF, '  SEE2  VERSION:  '.version); 

Writeln(OutF, '  INPUT  DATA  SET  NAME:  '.InFNam); 

Writeln(OutF,'  INPUT  DATA  SET  TITLE:  '.title); 

Writeln(OutF, ’  SetFNam-' .SetFNam, '  VFNam  -'.VFNam); 

Writeln(OutF, '  ACFNam**' .ACFNam, '  OutFNam*' .OutFNam) ; 

Writeln(OutF); 

Writeln(OutF, 'INITIAL  VALUES  OF  INPUT  VARIABLES:'); 

Out  pp ( pp , pp3, pp4 , ppe , conv , name , xndex .yndex ) ; 

Wrlteln(OutF) ; 

Writeln(OutF, 'ROW  VARIABLE:  ' ,name[ xndex] . '  varies  f rom' ,xmln: 9: 2, 

'  to' ,xmax:9: 1, '  in  steps  of ' ,xstep:9:2); 

WRITELN(OutF, 'FLOATING  VARIABLE:  ' .name [yndex] , '  (between  limits  of  , 
yll:7:2,'  and  ' ,yul: 9: 2, ' ) ' ) ; 

Wrlteln(OutF); 

Writeln(OutF, 'AIRCRAFT  TYPES:'); 

Wrlteln(OutF, '  ACID  TYPE  AX  AY  AZ'); 

For  1:*1  to  4  do 
begin 

If  (1<3)  then  iact:*Round(pp[ll])  else  If  (1*3)  then 
lact:>Round(pp4[ll])  else  lact;“Round(pp3[H ] ); 

Wrlteln(OutF, '  AC ,1 : l.acid [iact ] : 6, '  '  ,acname[iact] ,tabl [iact. 1] : 8: 1, 

tabl [iact ,2] :8: l.tabl [iact , 3] :8: 1); 
end ; 

Writeln(OutF) ; 

Writeln(OutF, '  -  OUTPUT  TABLE  - '); 

{ - START  COMPUTATIONAL  BLOCK - } 

lx:«0;  xval:“xmin-xstep;  yO: -pp [yndex] ; 

While  (xval<xmax)  do 

begin 

ix:-lx+l; 
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xval:"xval+x8tep;  ppAtxndex] :“xval*conv[xndex] ; 

WritelnC  >Computlng  for  ' ,name[xndex], '■* ,xval: 10; 1) ; 
Make_Con8l8tent(pp,pp3,pp4) ; 

{ - FIND  REFERENCE  PACQ - } 

pacqac [ 1 ] : -Fpacq ( pp , tab 1 ) ;  pacqac [ 2 ] : *pacqac [11; 
pacqref : “1 . 0-( 1 . 0-pacqac 1 1 ] )*( 1 . 0-pacqac [ 2 ] ) ; 

{-  -  ITERATE  PP4[YNDEX]  TO  OBTAIN  EQUIVALENCY  FOR  SECOND  PAIR  OF  AC - } 

nitt:"0;  n8top:**10; 

y ;-pp4[yndex] ;  ye:-y/conv[yndex] ;  dy:-10*dy8top*conv[yndex] ; 

While  (nltt<nstop)  and  (Ab8(dy)>dy8top*conv[yndex] )  do 
begin  {BEGIN  NITT  LOOP  } 
nltt;“nltt+l; 

IF  (KEYPRESSED)  THEN  BEGIN  WRITE(’CR  TO  CONTINUE’);  READLN(I);  END; 

IF  (debug>-2)  THEN  BEGIN 
WRITELN(OutF); 

WRITELN(OutF, '  -  BEGIN  ITERATION  NO.  •,NITT;3,'  -  ’); 

WRITELN(OutF, ’  VALUE  OF  ' ,nanie[yndex] , ’  IS  NOW’ ,y:10;4);  END; 
£old:>£new;  yold:*ynew;  ynew:-y; 

pacqac ( 3 ] ; -Fpacq(pp3, tabl ) ;  pacqac [4 ] ; *Fpacq (pp4, tabl ) ; 

IF  (debug>2)  then  WRITELN(OutF,  ’  FOR  ’,naine[yndex], '  “  *,ye:7:2,  '  kt', 

’  pacqac[4]-',pacqac[4]:9:5); 

£new;“1.0-( 1. 0-pacqac [3] )*(1.0-pacqac [4] ); 

If  (nltt>l)  then 
begin 

If  (ynew^old)  then  Writeln(Chr(7), ’ ! ICONVERGENCE  PROBLEM  ;  YNEW-YOLD-', 
YNEW:10;4,'  AT  Y-',Y:12;4) 

elae  If  (fnew-fold)  then 

Writeln(OutF,Chr(7), ' ! !  CONVERGENCE  PROBLEM:  fnew-fold-' ,fnew; 9: 5) 

else 

begin 

dfdy;“(fnew-fold)/(ynew-yold); 
dy : -(pacqref -f  new) /df dy 
end; 
end; 

y:*y+dy;  ye :-y/conv[yndex] ; 

If  (ye<yll)  then  ye:-yll;  If  (ye>yul)  then  ye;-yul; 
y : -ye*conv [yndex ] ;  pp4[yndex] ;-y; 

Make_Con8i8tent(pp,pp3,pp4) ; 

If  (debug>2)  then 
begin 

WRITE(OutF,'  NITT  YNEW  DY  FNEW  FOLD  DFDY'); 

WRITELN (OutF , ’  PACQAC [ 3 ]  PACQAC [ 4 ] ’ ) ; 

WRITE ( Out F, NITT: 4, YE; 10: 3, 'kt  ’ ,DY/conv[yndex] : 8; 3, 'kt ' ,FNEW; 10:4,’  ’ , 
FOLD;8:4,'  ',dfdy;9:4); 

WRITELN(OutF,pacqac[3] : 10:5,pacqac[4] ; 10:5); 
end; 

end;  {of  nitt  loop} 
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ENDITT; 

If  (nitt-nstop)  or  (y-yll)  or  (y-yul)  then 
begin 

Writeln(OutF) ; 

Writeln(OutF, ’EQUIVALENCY  SEARCH  FAILED  FOR  X  ,xval : 10: 4, ’  (  nitt.y- 
nitt:4,y: 12:4, ' 
end ; 

{ - end  computational  block - 

L2: 

If  (xval*xinin)  or  (debug>0)  then  Writeln(OutF, '  ' ,name[xndex] , ' 

name[yndex] , '  pacqref  pacq3  pacq4  ’); 

Writeln(OutF, '<<’ ,xval/conv[xndex] : 11: 3,y/convtyndex] : 10:2, 
pacqref : 9:4,pacqac[3] : 9: 4,pacqac [4] :9: 4, *  »' ); 
end;  {of  x  index  block} 

{pp[yndex] : “nomval [yndex] *conv [yndex] ;  } 

Writeln(OutF) ; 

WritelnCOutF, '*-*-*-*  END  OF  SEE2  ANALYSIS  *-*-*-*’); 

CloseCOutF); 

Writeln( 'Analysis  Complete  -  results  stored  in  ',0utFNam); 

Writeln( 'MENU: ' ) ;  WritelnC  1  “  print  results'); 

Wrlteln('  2  =  display  results');  WritelnC  CR  “  continue'); 

1:**0;  Readln(l); 

If  OO)  then 
begin 

If  (1“1)  then  AuxOutPtr :*LstOutPtr  else  AuxOutPtr:“ConOutPtr ; 

Assign(0utF,0utFNam) ;  Reset(OutF); 

il:-0; 

While  (not  Eof(OutF))  do 
begin 

Readln(0utF,s) ;  Wrlteln(Aux,s );  11:*'11+1; 

If  (il>“12)  then  begin  Delay(2000);  ii:“0;  end; 
end; 

Close(OutF) ; 


end ; 

Write('CR  to  continue');  Readln(l); 

1:=0;  WritelnC 'MENU  :'); 

WritelnC  1  -  NEW  ANALYSIS  (MODIFY  PARAMETERS)'); 

WritelnC  CR  -  QUIT'); 

Write('»');  Readln(i); 

If  (i”l)  then  goto  REDEF; 
end. 
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